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ABSTRACT 

Extreme ultraviolet (EUV, 13.6 eV < hu < 100 eV) and X-rays in the 0.1-2 keV band 
can heat the surfaces of disks around young, low mass stars to thousands of degrees 
and ionize species with ionization potentials greater than 13.6 eV. Shocks generated 
by protostellar winds can also heat and ionize the same species close to the star/disk 
system. These processes produce diagnostic lines (e.g., [Nell] 12.8 /xm and [01] 6300 
A) that we model as functions of key parameters such as EUV luminosity and spectral 
shape. X-ray luminosity and spectral shape, and wind mass loss rate and shock speed. 
Comparing our models with observations, we conclude that either internal shocks in the 
winds or X-rays incident on the disk surfaces often produce the observed [Nell] line, 
although there are cases where EUV may dominate. Shocks created by the oblique 
interaction of winds with disks are unlikely [Nell] sources because these shocks are too 
weak to ionize Ne. Even if [Nell] is mainly produced by X-rays or internal wind shocks, 
the neon observations typically place upper limits of < 10^^ on the EUV photon 
luminosity of these young low mass stars. The observed [01] 6300 A line has both a 
low velocity component (LVC) and a high velocity component. The latter likely arises 
in internal wind shocks. For the former we find that X-rays likely produce more [01] 
luminosity than either the EUV layer, the transition layer between the EUV and X-ray 
layer, or the shear layer where the protostellar wind shocks and entrains disk material 
in a radial flow across the surface of the disk. Our soft X-ray models produce [01] LVCs 
with luminosities up to 10~^ Lq, but may not be able to explain the most luminous 
LVCs. 



1. INTRODUCTION 

The photoevaporation of a protoplanetary disk by the extreme ultraviolet (EUV, 13.6 eV 
< hu < 100 eV) or Lyman continuum photons from the central star may significantly affect the 
formation and evolution of planets and planetesimals, and may be one of the important mechanisms 
for dispersing disks (Hollenbach et al 1994, Hollenbach, Yorke & Johnstone 2000, Clarke et al 
2001, Richhng, Hollenbach, & Yorke 2006, Alexander et al. 2006a,b, Alexander 2008a). EUV 
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photoevaporation occurs because the EUV photons create a 10^ K, ionized surface on the disk, and 
beyond about 1(M*/1 Mq) AU, where M* is the steUar mass, the thermal pressure of the gas is 
sufficient to drive a significant hydrodynamic flow out of the gravitational potential of the star and 
into interstellar space. 

Some of the most detailed models of the dispersal of disks around isolated low mass stars 
invoke viscous spreading and accretion on the inside ( < few AU) of the disk and EUV-induced 
photoevaporation on the outside (Clarke et al 2001, Matsuyama et al 2003, Ruden 2004, Alexander 
et al 2006 a,b). This combination has been invoked to explain gas-poor giants like Uranus and 
Neptune (Shu, Johnstone, & Hollenbach 1993), the rapid evolution of classical T Tauri stars to 
weak-lined T Tauri stars (Clarke et al 2001, Alexander et al 2006 a,b), the production of large inner 
holes such as exist in some sources (Alexander et al 2008b, Cieza et al 2008), and the migration 
and "parking" of giant planets (Matsuyama et al 2003, Lecar & Sasselov 2003, Veras & Armitage 
2004). 

X-rays from the star also significantly affect the disk. Glassgold et al (2004) show that hard 
X-rays can penetrate to moderate depths into the disk and produce sufficient ionization to maintain 
a vigorous magnetorotational instability (MRI, Balbus & Hawley 1991), at least in the upper layers 
of the disk (see also, Sano et al 2000, Stone & Pringle 2001). Chiang & Murray-Clay (2007) have 
recently expanded on this idea, using X-rays to stimulate MRI in the inner edge of a dusty disk, 
thereby eating away the disk from inside out. Alexander et al (2004a) have argued that X-rays 
by themselves do not lead to significant photoevaporation, but Gorti &: Hollenbach (2009) have 
shown that ~ 1 keV X-rays may increase FUV-induced photoevaporation rates by roughly a factor 
of 2, because X-ray ionization increases the electron abundance, which enhances the FUV grain 
photoelectric heating mechanism. More recently (Ercolano et al 2008, 2009, Gorti, Dullemond, 
&; Hollenbach 2009), it has become clear that a soft (0.1-0.5 keV) X-ray component can lead to 
significant photoevaporation rates. Glassgold, Najita, k. Igea (2007) and Meijerink et al (2008) 
have shown that X-rays partially ionize and heat the gas just below the EUV fully ionized layer, 
and that the X-ray heated gas achieves temperatures of order 1000 to 4000 K in a dense (n 10^ 
cm~'^) layer out to about 10-20 AU. Although they do not discuss X-ray photoevaporation, these 
temperatures, densities and radii suggest significant rates. 

EUV and X-ray photons around low mass stars, whose photospheres are too cool to produce 
a substantial EUV or X-ray flux, emanate from the accretion shock created by the impact of the 
accreting disk gas onto the stellar surface and/or from the hot plasma generated by magnetic 
activity on the stellar surface akin to (but much greater than) the Sun's chromosphere or corona. 
These two mechanisms heat plasma to temperatures >> 10^ K, and thereby produce significant 
EUV and X-ray luminosity. Alexander et al (2004b) argue that the EUV photons do not likely 
penetrate the accretion columns to irradiate the disk, and that, therefore, magnetic activity is a 
more attractive source for the EUV that shines on the disk surface. However, many accreting 
sources exhibit a soft X-ray component (e.g., Kastner et al. 2002, Stelzer & Schmitt 2004), which 
may arise from an accretion shock or is at least mediated by accretion flows (Preibisch 2007, Giidel 
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Sz Telleschi 2007). Soft X-rays are only somewhat more penetrating than EUV photons, raising the 
possibiUty that the geometry of the accretion streams (sometimes called "funnel flows") onto the 
star may also allow the escape of at least some of these hydrogen- ionizing photons. The hard ( ^ 1 
keV) X-rays likely arise from the magnetic activity (i.e., the chromosphere and corona). 

Whatever the source of EUV photons, they must still penetrate protostellar winds. Protostellar 
winds are thought to be driven by magnetohydro dynamic processes from the inner portions of 
accreting disks (e.g., Shu et al. 1994, Ouyed & Pudritz 1997). We show in §2 that these winds 
must have low mass loss rates, < 10~^ Mq yr"-*^, for EUV or soft ( < 0.2 keV) X-rays to penetrate 
them and to illuminate the disk surface beyond ^ 1 AU, where photoevaporation proceeds. Since 
accretion rates onto the central star are correlated with protostellar wind mass loss rates (Hartigan 
et al. 1995, White &: Hillenbrand 2004), this critical wind mass loss rate corresponds to an accretion 
rate of about 10"^ Mq yr"^ 

The main weakness in EUV photoevaporation models is the extreme uncertainty in the EUV 
photon luminosity ^euv of the central star. The EUV opacity of hydrogen is so high that a column 
of only ~ 10^^ hydrogen atoms/cm^ provides optical depth of order unity. Therefore, interstellar 
extinction prevents the direct observation of the EUV flux from young, low mass stars with disks. 
There are, however, observations of nearby, older, solar-mass stars, including the Sun, which provide 
a clue to the spectra from the far ultraviolet (FUV, 6eV < hu < 13.6 eV) to the X-ray of low mass 
stars due to their magnetic activity (Ribas et al 2005). These suggest that, very roughly, vF^ 
is constant for a given star from the FUV band through the EUV band to the keV X-ray band. 
Ercolano et al (2009) also discuss observations of flare stars which suggest magnetically-heated 
coronae on the stellar surfaces with a range of plasma temperatures resulting in roughly Fy oc 
power law EUV spectrum. Thus, one might estimate the magnetically-produced EUV luminosity 
of a low mass star by measuring either (or both) the 0.1-1.0 keV X ray luminosity or the 6 - 13.6 
eV FUV luminosity. For non-accreting but young (~ 1 Myr) low mass stars, the X ray and FUV 
luminosities tend to be of order ~ 10~^ L^o; (e.g., Flaccomio et al. 2003, Valenti et al. 2003), 
suggesting Lew ^ 10"'^ L0 or ^^euv ~ 10^^ EUV photons per second for a 1 M© star. 

Alexander et al (2005), based on earlier work of Brooks et al (2001), use FUV emission lines 
of various ions of the elements carbon, oxygen, nitrogen and silicon seen in T Tauri stars to try 
to estimate the distribution of emission measures as a function of T of the hot (~ 10^ — 10^ K) 
plasma. As the authors themselves point out, this method is fraught with difficulties, and, as a 
result, they can only constrain ^euv to range from 10^^ — 10^^ s~^ in young solar mass stars. 
The knowledge of the EUV luminosity is critical in predicting EUV-driven photoevaporation and 
determining whether it dominates disk evolution and explains the observed short (~ 1 — 3 Myr) 
lifetimes of disks around low mass stars. 

One way to measure ^euv is to observe emission lines produced by the heating and ioniza- 
tion caused by these photons on the disk surface. Such a measurement is important since ^euv 
determines the EUV photoevaporation rates, and therefore the EUV-induced dispersal times of the 
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gas and dust in these young, planet-forming disks. Given a disk illuminated by EUV photons, a 
tenuous, 10^ K, fully ionized surface is created by the photoionization of hydrogen. In effect, a sort 
of "blister HII region" is created above the bulk of the disk, which is mostly neutral molecular gas. 
Although this HII region contains very little mass (~ 10-'^ EUV /W^^ s-^]^/^Mj, where Mj IS a 
Jupiter mass), it can produce sufficient forbidden optical line emission (e.g., [SII] 6731 A and [Nil] 
6583 A, see Font et al 2004) or infrared fine structure emission (e.g., [Nell] 12.8 fim, this paper) to 
be observed. We note that [Nell] 12.8 fj,m is one of the strongest lines from HII regions associated 
with Giant Molecular Clouds, and, because neon is not depleted and its gas phase abundance rel- 
ative to hydrogen is quite well known, this fine structure line can also be used in these regions to 
measure or constrain the ionizing luminosity of the exciting star(s) (Ho & Keto 2007). 

There are two problems in using the emission lines from the HII surface to measure ^euv- 
Uncertainties in extinction, the gas temperature and the gas density make the optical lines a poor 
diagnostic of^EUV- The infrared fine structure lines are much better for this purpose, but they also 
can be produced by the heating and (partial) ionization of the neutral gas below the HII surface 
by penetrating X rays (Glassgold, Najita, & Igea 2007). In addition, they can be produced in high 
velocity (ionizing) shocks created by the protostellar wind. We discuss in this paper the relative 
contributions to the fine structure emission by the surface EUV heated layer, the subsurface X ray 
heated layer, and the wind shocks. However, even if its origin cannot be distinguished, the fine 
structure emission, e.g., [Nell] 12.8 iJ,m, gives a strict upper limit on ^euv- In addition, if arising 
from the EUV or X-ray layers, the [Nell] and other fine structure lines provide a measure of the 
density and temperature of the hot surface gas, and therefore directly probe some of the regions 
where photoevaporation originates (Alexander 2008). 

[Nell] 12.8 /im emission from young stars with optically thick disks was first detected using 
the high resolution mode of the IRS instrument on the Spitzer Space Telescope (Pascucci et al. 
2007, Lahuis et al. 2007, Ratzka et al. 2007, Espaillat et al. 2007), and is now found in over 
~ 50 sources (Giidel et al 2009). Some of these sources (~ 15) also show emission from the 
hydrogen recombination lines H(7-6)a and H(6-5)a, and only one source is detected in [NeIII]15/im. 
Observed line luminosities range from 10^"^ - 10^^ Lq. Follow -up, very high resolution ground- 
based observations of some bright [Nell] sources have resolved the line emission and observed line 
widths (~ 15 — 80km s~^), interpreted as emission arising from X-ray heated layers in Keplerian- 
rotating disks (Herczeg et al. 2007, Najita et al. 2009), EUV photoevaporative flows (Herczeg et 
al. 2007, Pascucci & Sterzik 2009) or from outflows associated with these sources (van Boekel et 
al. 2009, Najita et al. 2009). Correlations have been sought between the [Nell] luminosities and 
disk and stellar diagnostics such as X-ray luminosity (Pascucci et al. 2007, Giidel et al 2009) and 
mass accretion rates (Espaillat et al. 2007, Giidel et al 2009), but the data is inconclusive. The 
origin of the [Nell] emission, although widely attributed to disks, is still not definitive. 

This paper is motivated by the recent observations of [Nell] 12.8/im emission. We model disks 
illuminated by EUV and X rays, and present results for the infrared fine structure lines of Ar'^, 
Ar"'"'', Ne"*", Ne"'"^, N"*", N'''^, 0"*"+, S"*""*", and S"'"^"'', two infrared recombination lines of H"*", and 



- 5 - 



the optical forbidden line [01] 6300 A. We show that if the EUV layer dominates the emission, the 
infrared fine stucture lines diagnose ^euv and the shape (slope) of the EUV spectrum. We also 
show that measurements of [Nell] 12.8 /im and [Nelll] Ibfim are particularly good diagnostics of 
these parameters, being strong and relatively insensitive to extinction and changes in the plasma 
density n or temperature T. Our models of the X-ray layers, like the X-ray models of Glassgold 
et al (2007) and Meijerink et al (2008), produce [Nell] emission that, at least in some cases, is in 
accord with the observations. However, in a number of cases the X-ray heating mechanism seems 
insufficient to provide the emission (Espaillat et al. 2007, Giidel et al 2009), as we will also show 
in this paper. Shocks in the protostellar wind or an unseen EUV or soft (~ 0.1 — 0.3 keV) X-ray 
component may provide the origin of the [Nell] in these cases. Our models differ from Glassgold 
et al in that we treat the vertical structure of the disk consistently (that is, the gas temperature is 
not assumed to equal the dust temperature in calculating the vertical density structure), include 
EUV ionization and heating, include FUV photodissociation and heating, treat the X-ray heating 
somewhat differently, and include some additional significant cooling lines, such as [Nell] 12.8 
and [Aril] 7 fim. 

This paper complements earlier (Gorti & Hollenbach 2004, 2008) papers which examined the 
molecular and atomic fine structure emission from the neutral disk. In this older work, the fine 
structure lines treated focussed mainly on those with ionization potentials less than 13.6 eV, such 
as those of O, C, C+, S, Si, Si+, Fe, and Fe"*", although we did treat the X-ray ionization of 
some species in the predominantly neutral gas. In this paper we focus on species with ionization 
potentials greater than 13.6 eV, which are only found in the fully photoionized HII region surfaces 
of disks, in X-ray ionized, predominantly neutral gas, or in fast (^100 km s~^) shocks produced 
by the stellar wind. 

We organize the paper as follows. We discuss the restriction on the wind mass loss rate in 
order for the FUV, EUV and X-ray radiation from the protostar to penetrate the wind and shine 
on the disk surface in §2. Section 3 provides analytic estimates of the relation of the fine structure 
and hydrogen recombination line luminosities to ^euv^ the scaling of the emission from the X-ray 
layer to the X-ray luminosity of the central star, the [Nell] luminosity produced in wind shocks, 
and the [01] 6300 A luminosity possible from both the disk and from wind shocks. Section 4 shows 
the results of numerical models. Section 5 compares the results of our models to recent observations 
made by the Spitzer Space Telescope and several ground-based telescopes, and discusses the relative 
contributions of EUV, X rays and shocks to the observed [Nell], hydrogen recombination lines, and 
[01] emission. We conclude with a discussion and summary in §6. 

2. FUV, EUV AND X-RAY PENETRATION OF PROTOSTELLAR WINDS 

Although our protostellar wind model is influenced by the "X wind" models of Shu et al 
(1994), the main assumption we make is that the bulk of the wind mass loss rate originates 
from cylindrical radius to + /r^, where ~ 10^^ cm and f 1. Therefore, the model 
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applies also to other disk wind models (e.g., Ouyed & Pudritz 1997) where the bulk of the mass 
loss originates from the inner disk surface. We assume that / is sufficiently small that we can take 
rib as the average hydrogen nucleus density at the base of the wind without introducing significant 
error by assuming this constant density from to + fvy^. 

The mass loss in the wind, M^, arising from this geometry is given 

= 2[7r{rw + fr^)"^ - nrljubmnVw, (1) 

where mn = 2.3 x 10~^^ gm is the mass per hydrogen nucleus and is the wind velocity. The 
hydrogen nucleus column density N^, through the base of the wind, which the energetic photons 
must penetrate to reach the outer disk surface, is then given: 

Interstellar dust requires a hydrogen nucleus column of ^ 10^^ cm~^ to provide optical depth 
unity in the FUV. However, the dust lifted from the surface of the disk at the base of the wind 
is likely to have coagulated to much larger sizes than interstellar dust, and furthermore to have 
lower dust/gas mass ratios because of sublimation of the less refractory materials and settling of 
the refractory grains to the midplane (Dullemond & Dominik 2005). In fact, at radii of < 10^^ 
cm it is possible that all dust has sublimated. All these processes lower the dust cross section per 
hydrogen nucleus. Even if there is no dust (for example, if all the dust is sublimated), the gas 
provides FUV opacity and attenuates the FUV significantly for columns greater than about 10^^ 
cm~^. Assuming a minimum reduction in dust opacity relative to interstellar dust of a factor of 10, 
FUV will penetrate wind columns Ny^ < 10^^ cm~^. Dust also provides a source of X-ray opacity, 
which will be reduced from interstellar values by the effects of settling and coagulation. However, 
considerable opacity remains in the gas phase elements such as C, O and Ne. Gorti & Hollenbach 
(2004, 2008) estimate, using the cross-sections of Wilms, Allen & McCray(2000), that ~ 10^2 
cm~2 is required for 1 keV optical depth unity at disk surfaces. On the other hand, soft X-rays 
experience considerably more optical depth, and A'^ ~ 10^'' cm"^ provides optical depth unity for 
~ 0.2 keV X-rays. Therefore, in summary, M^, < 4 x 10~^ Mq yr~-^ is required for ~ 1 keV X rays 
to penetrate the protostellar wind, whereas soft X-rays can only penetrate when ^4 x 10"^*^ 
M0 yr~^. The penetration of the FUV likely occurs at mass loss rates considerably higher than 
4 X 10"^ Mq yr~^ because of dust sublimation and settling, but this number serves as a 
useful lower limit. 

A column Nyj of 10^" — 10^^ cm"^ of neutral hydrogen is totally opaque to EUV photons, since 
N[HI) ~ 10^'' — 10^^ cm~2 produces EUV optical depth unity. For EUV photons to penetrate 
the wind, the EUV photon flux Feuv must be sufficiently high to keep the base of the wind fully 
ionized, so that n{HI)/nb << 1 and N{HI) ^10^^ cm"^. This "Stromgren" condition 



Feuv > ar,Hnlfry,, 



(3) 
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can be rewritten as 

A> in-io f^^ ^^^V '^Euv A^^V \ -1 

Muj<8 X 10 — — Tj r T yr \ (4) 

In other words, the mass loss rate has to be less than the RHS of this equation for EUV to penetrate 
the base of the wind and illuminate the outer disk surface beyond 1 AU. 

Wind mass loss rates are hard to measure "directly" from the observed optical line emission 
(e.g., [SII]) seen in their jets. The derived mass loss rates from these optical lines depend on knowing 
the gas temperature, the gas density, and the ionized fraction - all of which are quite uncertain. 
"Indirect" methods rely on measuring the momentum in swept up circumstellar gas. This method 
is also approximate since it requires an estimate of the wind speed, the duration of the mass loss 
episode, and the conversion factor of CO luminosity to mass. Likewise there are uncertainties in 
observationally determining the mass accretion rate Mace onto the central star. These uncertainties 
create a spread in the constant k of proportionality, but it is generally agreed that the wind mass 
loss rate scales with the mass accretion rate, ~ kMacc- The constant k has been estimated 
from ~ 0.01 (Hartigan et al 1995) to ~ 0.1 (White & Hillenbrand 2004). The Shu et al (1994) X 
wind model predicts values somewhat higher than 0.1. White and Hillenbrand point out that there 
seems to be considerable intrinsic scatter in the ratio of wind mass loss rate to mass accretion rate 
from source to source. 

Roughly then, if we take ~ 0.1 Mace, the FUV and ^ 1 keV X-rays penetrate the wind 
when Mace ^4 x 10~^ M.q yr~^, whereas the EUV and soft (~ 0.2 keV) X-rays penetrate the wind 
when the accretion rate has dropped to Mace ^ 8 x 10~^ Mq yr~^. Hartmann et al (1998) show the 
evolution of Mace for young, solar mass stars. With order of magnitude dispersion. Mace is roughly 
10~^ M0 yr~^ at 1 Myr, and drops rapidly on Myr timescales. Thus, FUV and 1 keV X rays 
illuminate the disk surface nearly as soon as the epoch of heavy accretion of material onto disk and 
star from the natal cloud core has ceased. However, EUV and soft X-rays may not illuminate the 
disk surface until roughly 1-2 Myr has elapsed from that time. 

If one wishes to observe a disk whose ionized fine structure lines are not produced by EUV 
and soft X-rays, one should select sources with Maec ^8 x 10~^ yr~^. If [Nell], for example, is 
detected in sources with 8 x 10~^ M0 yr"-"^ < Maec ^ 4 x 10"'^ yr~^, then hard (~ 1 keV) X-rays 
or possibly wind shocks may be implicated. If [Nell] is detected in sources with Mace > 4x10-^ Mo 
yr~^, then protostellar wind shocks almost certainly provide the origin. If one wishes to observe 
sources illuminated by both X-rays and EUV, and therefore containing, for example, [Nell] emission 
from an EUV-produced HII surface and also from an X-ray-produced partially ionized deeper layer, 
then one should observe sources with Maec ^8 X 10-9 Mq yr-^ Interestingly, the [Nell] sources 
have been detected in [Nell] emission in this entire range of Maec, suggesting a wide range of origin 
of the [Nell] [see Espaillat et al (2007) and Giidel et al (2009)]. 
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3. ANALYTIC MODELS OF EMISSION LINES DIAGNOSTIC OF EUV AND 



There are basically three types of lines diagnostic of EUV and X-rays incident on disks or of 
fast, ionizing shocks in protostellar winds: hydrogen recombination lines, optical forbidden lines like 
[01] 6300 A, [SII] 6713 and 6731 A, and [Nil] 6583 A, and the infrared fine structure lines of ionized 
species whose ionization requires photons more energetic than 13.6 eV. We discuss here the last 
two, and leave the discussion of hydrogen recombination lines for §3.2. With the possible exception 
of [01] 6300 A, the optical lines likely arise in the completely ionized HII region at the surface of 
the disk or in fast, ionizing shocks because the lines typically lie > 20, 000 K above ground and 
are excited mainly by electron collisions; these regions have higher temperatures (~ 10"^ K versus 
several thousand K in the X-ray heated region) and generally higher electron densities than the 
X-ray layers. In addition, most of the optical lines arise from ionized species whose abundances 
peak in the completely ionized HII gas as opposed to the mostly neutral X-ray layers (the notable 
exception being [01]). The infrared fine structure lines from high IP (> 13.6 eV) species typically 
lie ~ 300 — 1000 K above ground and therefore are not sensitive to temperature for temperatures 
above about 300 — 1000 K. These lines may come from either the EUV- heated HII region, the X-ray 
heated region, or shocked regions and we show below that the relative EUV versus X-ray photon 
luminosity from the central star determines which of these two regions will dominate the emission. 
We focus in this section on the infrared lines, because Font et al (2004) and Meijerink et al (2008) 
have discussed the optical emission from the EUV and X-ray heated layers. However, we do include 
the [01] 6300 A line in our analysis, because our [01] luminosities from the X-ray layer differ from 
the Meijerink et al values, and because other researchers have not been able to match the observed 
luminosities in this line (Hartigan et al 1995). We also include a discussion of the infrared hydrogen 
recombination lines that have been observed from these star /disk systems. 



Consider an axisymmetric disk described by cylindrical coordinates r, z. If /euv is the fraction 
of ionized photons from the star absorbed by the disk, then the Stromgren condition is: 



where the electron density Ue is a function of r and z but is negligible below zjf, the ionization 
front, and where we have ignored dust attenuation in the HII surface region above zip. We will 
justify the neglect of dust post facto below, as well as show that /euv ~ 0.7. In Eq. (5) and 
To are the inner and outer radii of the disk, and ar,H is the case B recombination coefficient for 
electrons with protons {ar^n = 2.53 x 10~^^ cm'^ s~^ at T = 10^ K, Storey & Hummer 1995). 

For simplicity we treat the specific example of a simple two level fine structure system like 
[Nell] 12.8 /im. Then, for a transition t we have that the escaping line luminosity Lt from the disk 



X-RAYS INCIDENT ON DISKS OR OF SHOCKS 



3.1. Fine Structure Lines from the HII Surface (EUV Layer) 




(5) 
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IS given: 




(6) 



where 74 is the colHsional excitation rate coefficient of the transition, AE'j is the photon energy, 
Uecr^t is the critical electron density for the transition, and n{i) is the density of the ionized species 
i that produces the transition. Note that here we have included the fact that half the emitted 
infrared photons are directed toward the disk midplane, where they are absorbed by the (assumed) 
optically thick disk, and half escape. If we set n{i) = Xsf{i)ne, where Xg is the abundance of species 
s (in all ionization states) in the EUV layer, and f{i) is the fraction of that species in ionization 
state i, then we can write using Eq.(5): 



We note that in taking Xsf{i) out of the integral in Eq.(6) we implicitly assume that in Eq.(7) 
Xsf{i) is the density weighted average of this product in the EUV layer. Thus, we find that if 
Tie < Uecr^t and if f{i) does not depend on ^euv (e-g, for the dominant ionization state where 
/(i) ~ 1), the line luminosity Lt is directly proportional to the ionizing photon luminosity ^euv 
of the central star. If EUV is the sole excitation source, the measurement of Lt directly measures 
the uncertain parameter ^euv if we have a knowledge of Xs and /(«). The main unknown is the 
fraction f{i) in a particular ionization state i, since the gas phase abundance of an element can 
often be estimated from observations of HII regions. For example, neon often is found as Ne"*" and 
Ne"*"*". We therefore rewrite Eq. (7) as 



where Ct = ^t^Et^s / {"^otr^H) ■ Note that Ct and necr,t are known quantities and therefore constants 
in the equation. If the transition is from the dominant ionization state of the species, then /(i) ~ 1. 
Our modeling of disks suggests Jeuv ~ 0-7- Ct is therefore the main constant of proportionality 
that links the observed IR line luminosity to the EUV photon luminosity; for rig < necr,t and taking 
the dominant ionization state, Ct is the energy per absorbed EUV photon that emerges from 
the disk in the fine structure transition. Table 1 lists the value of Ct for the various transitions 
considered in this paper, along with necr,t and the assumed Xs that is used in Ct- 

The above equations show that it is important to estimate the electron density in the ionized 
surface of the disk at the radius where most of the emission in a given line is produced. We follow 
the results of Hollenbach et al (1994). Because the HII region surface is isothermal, the electron 
density at a given r decreases from zje upwards as 





(8) 




(9) 
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where H is the isothermal scale height of the 10^ K gas given by 



for r < r„, and where r„ is given by 



3/2 

H{r) = rAyJ (10) 



^5 ^ 7 ( ) AU. (11) 



1 M© 



Note that in Eq.(9) both zip and H are functions of r. The gravitational radius Vg is where the 
sound or hydrogen thermal speed is equal to the escape speed from the gravitational potential of the 
star. For r > Vg, where the gas is freely evaporating, the effective height of the disk is ~ r, since 
the initially vertically flowing gas turns over on radial streamlines by the time z = r. Hollenbach 
et al also showed that the base electron density at zjp falls off radially as a power law. 

n.(^^,.)-5xlO^(^)"'(rs)"'^(i^)''' c,.- fo„<., (12) 

and 

/r„\ 5/2 

ne{r,ziF) '^ne{rg,ziF) [ — ] for r > r^. (13) 



Note that if ne{r,zjF) < Uecr everywhere, then the amount of luminosity L{r) from a logarithmic 
interval of r is proportional to the volume emission measure at that r, or to n^r^H. For r < Vg 
we then obtain L{r) oc r^/^ whereas for r > we obtain L{r) oc r~^. In other words, the line 
luminosity originates mostly from r ^ rg. If the electron densities in the very inner regions exceed 
Uecr but are less than riecr at r^, the luminosity is relatively unaffected, and our conclusion does not 
change. However, if ne{rg, zjp) > Ugcr, then the line luminosity will drop. The line emissivity will 
be suppressed by a factor of approximately ne{rg, zip) /near at Vg. However, now L(r) oc rieV as 
long as He > Uecr, SO that L{r) oc r^/^ beyond Vg until the density drops below the critical density, 
when it reverts to its former dependence. The luminosity of a low critical density species will 
therefore originate from the radius 

neirg^Y^^^ (14) 

where n^irmax^ziF) = Uecr- As a result, 

L oc ne{rmax,ziF)r%,^^ oc ^%ly (15) 

as long as ne{rg, zif) > n^cr and f{i) is constant. 

As examples, consider [Nell] 12.8 fim and [SHI] 19 fim. The critical density for [Nell] is given 
f^ecr^Neii] = ^2i/72i ~ 6 X 10^ cm~'^, where we have taken A21 = 0.00859 and a collisional 
de-excitation rate coefficient at 10^ K of 721 = 1.355 x 10~^ cm^ s~^ (Griffin et al. 2001). Thus, 
comparing this critical density with the electron density at rg [Eq. (12)], we see that [Nell] is 
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subthermal at Vg typically, and that therefore it mostly originates from rg and tracks $ euv linearly 
as long as f{Ne~^) is constant. However, [SIII] has a critical density of approximately 5 x 10^ cm~'^ 
which is about 10 times less than ne{rg, Z[f) for a solar mass star with ^euv ~ 10"^^ s~^(see Eq. 
12). Therefore, 

'''max — W'^^^Tg — l.brg. The [SIII] line luminosity is down by a factor 10^^^ ~ 6.3 
from the value it would have had if it had been subthermal at r^, rather than the factor of 10 
drop at Tg, because most of the emission comes from further out where there is more volume. 
The luminosity in the line will not linearly track ^euv because of the significant (and variable 
with ^Euv) collisional de-excitation of the transition. In fact, as long as ne{rg, zjf) > necr end 
f{S~^^) is not dependent on ^euVj the luminosity in the line scales as ^euv ^ shown in Eq. (15). 
However, we find in our numerical analysis that f(S~^^) does depend significantly on ^euV: and 

3/5 

the ^EUV dependence is not seen (as we show below in our numerical models). 

We have shown above that the fine structure emission from the HII surface region arises from 
~ Vg as long as ne{rg, zjf) < n^cr- However, this conclusion and the important analytic derivation 
of the line luminosities [Eq. (8)] both require that the surface HII region be ionization-bounded 
(i.e., that there be a neutral layer underneath the completely ionized surface). The minimum mass 
of gas required to fully absorb the incident EUV photons and create an ionized surface region with 
a neutral midplane region for r < r^, is given: 

Mi{jj{min) = I 2dz I mHn(,{r, z)27rrdr, (16) 

J zjp J 

where mH is the mass of the ionized gas per electron ('--^ 2 x 10~^^ gm). The z integral can be 

/ X 3/2 

approximated as 2H{r) = 2rg ( f- j with the electron density fixed at the density at the ionization 

front ne{r,ziF)- The density ne{r,zjF) falls as r~^/^ for r < rg. Therefore, for r < Vg the mass is 
mostly at rg. Performing the integral, we obtain 

MHiiimin) ^ 2nrlmHne{Tg, zif) ^ lO^^o (^^^^^ j [j^j Mq. (17) 

Thus, assuming ne{rg) < necr,t, an extremely small gas mass, of order 10"'' Jupiter masses, inside 
of rg will provide the luminosities given by Eq. (8), using the values of Ct given in Table 1. These 
lines, then, are very sensitive diagnostics of the presence of trace amounts of gas at radii of order 
1-10 AU in disks. If there is less mass than MHii{min) at rg, then the resulting luminosities will 
be reduced by a factor M h 1 1 / M h i limin) . 

We now address whether dust extinction is important in the HII surface region. Since H{r) oc 
r^/^, the 10^ K surface of the disk is fiared. Most of the emission comes from rg or beyond as shown 
above. Using Eqs. (11 and 12), the attenuating column at rg is roughly 



Natt - n^ivg, zjf) xrg^5x 10^^ [-^^) (18) 



The dust at the surface of the disk is expected to have less opacity than interstellar dust, so that 
the column for EUV optical depth unity is greater than 10^^ cm~^. The above equation shows that 
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dust will not be important until ^euv > 10 s~ for a 1 Mq star. Low mass stars do not produce 
such high EUV luminosities. Therefore, our neglect of dust opacity in the HII surface is justified 
in disks around low mass stars. 

Finally, we estimate the fraction /euv of the ionizing photons emitted by the star that are 
absorbed by the disk. Most of the absorption (and consequent emission, as shown above) comes 
from r ~ or beyond, and here H ^ r. As we will show below in our numerical disk models, 
the underlying neutral disk also has considerable height. At r = 10 AU, zjp — 7.5 AU (see §4.2). 
Thus, the disk is opaque to EUV photons from the midplane to an angle from midplane of about 
40 degrees. In addition, the recombining hydrogen in the ionized gas above 40 degrees also absorbs 
EUV photons. Therefore, we estimate that the fraction Jeuv of EUV photons absorbed by the 
disk is about 0.7. A detailed hydro dynamical study is needed to more accurately determine this 
fraction. 

As our prime examples for this analytic analysis, consider specifically the case of the [Nell] 
12.8 /xm, [Nelll] 15 fim, and [Aril] 7 fiui luminosity emerging from a young disk. We choose these 
lines because [Nell] and [Nelll] have been observed, and [Aril] is predicted to be the brightest 
of the unobserved lines (see Table 1). In addition, they all have high critical densities so that 
ne{rg,zjF) < fleer for these lines as long as ^euv < 10^^ — 10^^ s~^. Assuming this condition is 
met and substituting Jeuv = 0.7 and the atomic constants into Eq. (8), we obtain 



Recall that f{Ne'^) is the fraction of neon in the singly ionized state in the region near rg which 
produces most of the emission. Luminosities greater than about 10"'' Lq are detectable from 
nearby (< 100 pc) sources by the Spitzer Space Telescope, as long as the line to continuum ratio is 
sufficiently large to enable detection of the line above the bright continuum. 

The effect on IR luminosity caused by holes in disks. The above analytic results apply for a 
disk that extends inwards to r < ~ 7 AU from the star. However, disks have been observed 
with inner holes, devoid of dust, that extend to > rg (e.g. Najita et al. 2007, Salyk et al 
2009). Regardless of the cause of these holes, if gas is absent inside of rj and > Vg, the disk 
vigorously photoevaporates at r^, a process which evaporates the disk from inside out (Alexander 
et al 2006a, b). Alexander et al show that the flux of EUV photons striking the inner wall of the 
disk creates a thin (thickness ~ H, the scale height of the neutral disk at rj) ionized layer. The 
Stromgren condition gives the electron density in the layer: 




(19) 




(20) 




(21) 



Few 




(22) 
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Assuming H ~ O.lrj (Alexander et al. 2006a, b), we obtain 

ne(r,) 2.5 X 10^ (^^^^j (^^— j cm'^. (23) 

Note that the electron density decreases as the inner hole size increases. If rj >> Vg, this leads to an 
increase in the luminosity of low critical density lines with respect to high critical density lines [see 
Eq. (8)]. However, for lines whose critical densities are larger than the electron density at rg, the 
presence of a hole of size >> rg does not appreciably affect the IR line luminosity. Essentially, 
this arises because the IR line luminosity is proportional to the number of EUV photons absorbed, 
fEUV^EUVi and this remains constant, with fEUV ^ 0.7, regardless of r,. Therefore, the IR line 
luminosity tracks ^euv as presented in Eq. (8). 



3.2. Infrared Hydrogen Recombination Lines from the HIT Surface 

The infrared hydrogen recombination lines can be analytically determined by noting that the 
luminosity in a given line produced by the transition n„ to rii is given 

POO l-To 

Lui = ttuiAEui / dz 27rrnldr, (24) 

J Zip Jri 

where aui is the rate coefficient for recombinations through the — ni transition, and AEui is 
the energy of the photon produced in this transition. Clearly, the hydrogen recombination lines 
also count EUV photons (see Eq. 5) and could be used to measure ^euv- However, hydrogen 
recombination produces weak IR lines compared to the fine structure lines such as [Nell] if the 
electron density Ue is less than the critical density Uecr of the fine structure transition, as can be 
seen by taking the ratio of the predicted line luminosities 

Lui o.uiAEui , . 

(25) 



L{NeII) -flNeII]AENeIlXNef{Ne+) 

The hydrogen recombination lines we are most interested in are the 7-6 (Humphreys a) and 6-5 
(Pfund a) at wavelengths of 12.37 and 7.46 fim respectively. These two lines have been reported 
observed toward stars with disks (Pascucci et al. 2007, Ratzka et al. 2007). Substituting the atomic 
constants for these transitions we obtain predicted ratios for the EUV-induced surface HII layer: 

^'^^ 0.008 (26) 



L{NeII) 

and 



0.02 (27) 



L{NeII) 

The observed ratios are close to unity! The predicted ratios are small because of the low ratio of the 
radiative recombination rate aui of hydrogen to the electronic collisional excitation rate of [Nell] 
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l[Neii]- Thus, we predict that these IR hydrogen recombination hnes must originate from another 
source if the [Nell] originates from the HII region surface of the disk. One place which would 
provide copious recombination line emission without producing even more [Nell] emission would 
be very high electron density regions, Ue » "^ecr, [Afe//]- these regions the [Nell] is suppressed 
relative to the hydrogen recombination lines due to the collisional de-excitation of the upper level 
of the [Nell] transition. Therefore, these observed recombination lines are possibly produced in 
very dense plasma very close to the star, in the stellar chromosphere, the accretion shock, or an 
internal wind shock if it is both high speed {vs ^ 100 km so that it produces ionized hydrogen) 
and occurs so close to the wind origin ( < 1 AU) that the postshock density is high enough to 
suppress [Nell] relative to the recombination lines. In any of these cases, the prediction is that the 
H recombination lines will be much broader (> 100 km s~^) than the [Nell] lines (~ 20 km s~^) in 
face-on disks. 



3.3. IR Fine Structure Lines from the X-ray Heated and Ionized Subsurface Layer 

Glassgold et al (2007) and Meijerink et al (2008) have presented models of the [Nell] 12.8 
fim emission and emission from other lines, such as [Nelll], [01], [SIII] and [SIV], produced in the 
X-ray heated layer that lies just below the ionization front created by EUV photons. This layer is 
predominately neutral, Xe ~ 0.001 — 0.1, depending on r and z, but with a higher ratio of Ne'''/Ne. 
Typically, the [Nell] emitting layer has T ~ 1000 — 4000 K. We also in §4 present numerical 
results from our models of the X-ray-induced fine structure emission, and in §5 discuss differences 
between our models and those of these authors. Here, we present a simple analytic estimate 
of the strengths of the fine structure transitions in X-ray-illuminated regions. These estimates 
are more approximate than those presented above for the EUV-dominated regions because of the 
uncertainties in estimating the gas temperature in this mostly neutral gas illuminated by a spectrum 
of X-ray photons. Nevertheless, they provide insight into the X-ray process, and into the relative 
strengths of X-ray induced fine structure emission in the X-ray layer as opposed to that produced 
by EUV photons in the surface EUV layer. 

The simplest derivation arises if we assume our "hard" X-ray spectrum, dominated by 1-2 keV 
photons. These photons are sufficiently energetic to ionize the K shell of Ne, and the ionization of 
Ne is dominated by direct X-ray photoabsorption, and not by collisions with secondary electrons. If 
we make the assumption in the X-ray layer that the atomic Ne absorbs a fraction f^^ of all ~ 1 keV 
X-rays, and that Ne"*" radiatively recombines with electrons, and we assume that rie < ?^ecr,[7Ve//]) 
we obtain in a manner completely analogous to the EUV layer's Eq. (5-7): 

,X _ fNeII]^E[NeII]f^efx^X 

^^^''^ ~ • ^^^^ 

Here, Tj^g//] is the collisional rate coefficient for [Nell] by electrons in the X-ray layer (that is, it is 
only different from 'y[Neii] in Eq. (7) because the X-ray layer is cooler than the EUV layer), is 
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the ~ 1 keV X-ray photon luminosity of the central source, fx is the fraction of X-rays absorbed 
by the disk in the X-ray layer, and ar^Ne is the recombination rate coefficient of Ne"*" with electrons 
in the X-ray layer. The fraction of ~ 1 keV photons absorbed by neon, f^^, is approximately the 
neon cross section at 1 keV divided by the total X-ray absorption cross section at 1 keV; using 
Wilms, Allen, & McCray (2000), we obtain f^^ ~ 0.21. We see that L-^^^jj^ scales linearly with 
^x, just as L[]\j^ijj scales linearly with ^euv in the EUV layer. The ratio of the [Nell] luminosity 
from the EUV layer to that in the X-ray layer is then given 

'-[Nell] yilkell]) \InJx) \ "r.H ) V, 7 

where f{Ne~^) is the fraction of neon that is singly ionized in the EUV layer. We take T ~ 10^ K 
for the EUV layer and Tx ~ 2000 K for the X-ray layer to estimate the recombination coefficients, 
and assuming that fx is approximately the height of the layer which becomes optically thin to 1 
keV X-rays from the star (roughly ~ 10^^ cm~^, or a column of about 10^^ to the star) divided 
by r, or fx ~ 0.25 (see Fig. 5). The EUV layer is more flared, hence fsuv ~ 0.7. In the EUV 
layer f{Ne^) ~ 1. Substituting into Eq. (29), we obtain 

"^[^^^^1 - 2 X 10"3e^ f ^^'j . (30) 



^[Nell] \ J 

The 1 keV X-ray photon luminosity ^x from a typical source is of order 10^^ photons s~^. The 
EUV luminosity ^euv is generally of order 10^^ photons s~^. Therefore, assuming the 1 keV X- 
rays are absorbed in regions with Tx > 1000 K, the [Nell] luminosity is expected to be marginally 
dominated by emission from the X-ray layer as opposed to the EUV layer. We show below in our 
numerical work that L^j^^jjj/ L^^j^ ~ 0.6 when ^euv = 10^^ s~^ and ^x — 10^^ s~^ and when 
the EUV spectrum is such to produce more Ne"*" than Ne^"*" in the EUV layer; this result agrees 
with Eq. (30). Note that ^euv = 10^^ s~^ and ^x ~ 10^^ s~^ corresponds to Leuv — Lx- In 
other words, if the central star emits the same EUV and X-ray luminosity, and the EUV has a 
soft spectrum which produces more [Nell] than [Nelll], there will be roughly 2 times more [Nell] 
luminosity arising from the X-ray layer than from the EUV layer. As we will show in §4, where we 
present results from our detailed numerical models, this conclusion that X-rays are more efficient at 
producing [Nell] emission does not depend strongly on the X-ray spectrum for reasonable choices 
of the spectrum. If we adopt a softer spectrum, the ionization of Ne is dominated by secondary 
electrons because most of the X-rays are absorbed by He, C or O. The gas is also hotter because 
there is more heating per unit volume due to the higher cross sections for softer X-rays. The net 
effect is that the [Nell] luminosity does not change much for fixed X-ray luminosity, even as we vary 
the X-ray spectrum. In the case of the "hard" X-ray spectrum, the reason X-rays are somewhat 
more dominant than the EUV is because for high temperature (T > 1000 K) gas such as exist in 
both the HII region and the X-ray heated region, the luminosity in the line depends mainly on 
the number of Ne"*" ions times the electron density. In the HII region, the vast number of EUV 
photons are used ionizing H, an extremely small fraction of the photons are used ionizing Ne, and 
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therefore the number of Ne"*" ions times the electron density is small (the xi\[ef{Ne~^) factor in Eq. 
(29); however, we explicitly include the relatively large fraction, f^^ = 0.21 of 1 keV X-ray photons 
that directly ionize Ne and lead to a large product of electron density times Ne"^ ions in Eq. (29). 
Therefore, no X]sfef{Ne~^) term appears in the denominator. Another way of understanding this 
result is that, although the EUV layer is completely ionized with f[Ne^) ~ 1 and Xg ~ 1, the 
EUV layer has much less column of Ne^ because H and He rapidly absorb the EUV photons; the 
penetrating X-rays partially ionize a much larger column. 



3.4. Shock Origin of Ionized Infrared Fine Structure Lines 

Lahuis et al (2007) and van Boekel et al (2009) discuss the possible origin of [Nell] emission 
from shocks generated by protostellar outflows. In order for a shock to produce significant [Nell] 
emission, the shock must ionize most of the preshock gas in order to produce high quantities of Ne^ 
and electrons. The fraction of preshock gas that is ionized by the shock is a very sensitive function 
of the shock speed (e.g, Hollenbach & McKee 1989, hereafter HM89). HM89 show that Vs Z 100 
km s""*^ is required to ionize most of the H and Ne and that the [Nell] emission rises very sharply 
with Vs and then plateaus above f s ^ 100 km s~^. This suggests that any possible shock must 
originate from internal shocks in the protostellar wind (which has terminal speeds of ~ 200 km 
s~^) or from the protostellar wind overtaking much slower moving outflow material. It is unlikely 
to originate from the shock produced by the wind striking the disk, since this shock is so oblique 
that the (normal) shock speeds are typically < 20 km (Matsuyama, Johnstone & Hollenbach 
2009). The total emission per unit area Ft from the postshock region of a radiative shock is given: 

1 „ 

Ft = -mHnoVs, (31) 

where mn is the mass per hydrogen nucleus and no is the hydrogen nucleus number density of the 
gas. The numerical results of HM89 can be approximated for the emission per unit area of [Nell] 
for shocks with Vs ^ 100 km s^ : 

^INell] - ( -I , no ) ^t, (32) 

where the dependence on density at high density arises because of collisional de-excitation of the 
upper state of [Nell] in the postshock gas. Assume that a fraction fsh of the protostellar wind 
shocks at speeds Vg 100 km s~^ with preshock density uq. It follows that the [Nell] shock 
luminosity is 

LfLin = Ifs^^^-s ( ) - f ) fs^M.s Le, (33) 

^ V ^ + lO-i cm-;^ / + 10-1 cm-^ / 

where M-s = Mw/W~^ Mq yr~^ and f > 100 km s~^. Therefore, if the protostellar wind mass 
loss rate M^, > 10~^ Mq yr~-^, Vg ^ 100 km s~^, fsh'^^, and no ^ 10^ cm~^, the [Nell] luminosity 
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produced in internal wind shocks may be comparable to or greater than the luminosity produced 
in the EUV or X-ray layer of the disk. 

Van Boekel et al (2009) argue this may be the case in T Tau S. Here, the observed [Nell] 
luminosity is i[Are//] ~ 10~^ L0. Prom the above, this would require, for example, fsh ~ !> 
2.5 X 10~^ Mq yr"^, and no ^ 10'* cm~^. The preshock density can be estimated with 
knowledge of the distance of the shock from the star, r^fi, and My^. The preshock density (the 
density in the wind at rsh) is given 

no = 2.5 X W^M.sriifa^ cm-^ (34) 

where ris = r^/j/lO^^ cm and fa is the fraction of Att sterradians into which the protostellar 
wind is collimated. Van Boekel et al (2009) measure an extent of the emission from T Tau S of 
approximately 160 AU. If we assume Vsh = 160 AU and = 1, presumably upper limits, we obtain 
a lower limit to 71q ^ 1.5 x 10^ cm~^. Note that for > 2.5 x 10~^ yr~^ and for our specific 
assumptions on rgh and /n , the preshock density is so high that the [Nell] luminosity is independent 
of M^; the luminosity from the shock saturates once the emitting [Nell] in the postshock gas is 
in LTE. Therefore, although it pushes parameters a bit uncomfortably, if T Tau S has protostellar 
mass loss rates > 2.5 x 10"'' Mq yr~^, it is possible that internal wind shocks produce the observed 
[Nell] luminosity. Note that the "dynamical time", Vs/rsh ~ 10 yrs, is marginally consistent with 
the observations of no significant time dependence since 1998 (van Boekel et al 2009). 

More recently, Giidel et al (2009) have assembled [Nell] data from a large number of sources 
and have plotted L[7Ve//] versus Mace- For low values of Mace ^ 10~^ Mq yr~^ the [Nell] luminosity 
is nearly independent of Mace, and is typically of order 3 x 10~^ Lq. However, for higher mass 
accretion rates, and in particular, for all the sources with known outflows or jets, L[Are//] increases 
with increasing Mace (arguably linearly, but with a lot of scatter). We present these observational 
results from Giidel et al and compare them to Eq. (33) in §4.4. The correlation of -Z>[jve//] with 
Mace suggests either that the higher luminosity {L^^^^j^ ~ 10~^ — 10~^ Lq) sources may originate 
in protostellar shocks or from EUV or soft X-rays produced by the accretion of disk material onto 
the star. However, in the latter case, these photons must penetrate the disk wind, which seems 
unlikely. 

3.5. [OI] 6300 A Emission from Young Stars with Disks 

[01] 6300 A emission is often observed in young low-mass stars with disks and outflows (e.g., 
Hartigan et al 1995). Two velocity components are seen: a high velocity component "HVC" and 
a low velocity component "LVC". Hartigan et al argue that the HVC comes from shocks in the 
protostellar wind, similar to our discussion above of internal shocks. The typical velocity of this 
component is ~ 100 - 200 km s-^ and the [01] 6300 A luminosity is ~ 10"^ - 10"^ Lq. HM89 
show that for Vs ~ 100 km s~^ the [01] 6300 A emission from the shock is about 10 times more 
luminous than the [Nell] 12.8 //m emission for no ^ 10^ cm~^, with even higher ratios at no > 10^ 
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cm ^ because [01] 6300 A does not coUisionally de-excite as readily as [Nell] 12.8 fim. Therefore, 
in agreement with Hartigan et al, we find that mass outflow rates of > 10^^ yr^^ can produce 
the most luminous [01] 6300 A HVC sources (see Eq. 33). 

Hartigan et al (1995) attribute the LVC to [01] emission emanating from the disk surface, 
probably in a relatively slow outflow since the emission is observed to be slightly blueshifted (~ — 5 
km with great dispersion). However, there are red and blue wings extending to ±60 km s~^ 
in the LVC, presumably due to a combination of Keplerian rotation and outflow. The [01] 6300 
A luminosity in the LVC ranges from ~ 10-^ to ~ 10-3 L©, with "typical" values ~ 10"^ Lq 
(Hartigan et al 1995). The exact origin of the LVC [01] emission, and its heating source, remains 
a mystery. In Section 4 below we present our numerical model results for [01] emission from the 
EUV layer and the X-ray layer. In agreement with previous work by Font et al (2004), we find 
that the EUV layer can only provide [01] 6300 A luminosities ^10"^ Lq. Meijerink et al (2008) 
were able to produce [01] luminosities as high as 5 x 10~^ Lq in their models of the X-ray layer. 
Therefore, they found it very difficult to explain the most luminous [01] LVC sources, but were 
able to produce luminosities in accordance with many of the observations. However, our more 
detailed models with a similar X-ray spectrum as that assumed by Meijerink et al produce lower 
[01] luminosities, primarily because we obtain lower gas temperatures in the X-ray heated layer 
(see below). However, if we use a softer X-ray spectrum such as the one proposed by Ercolano et 
al (2009), we do obtain luminosities of order 10~^ Lq. In summary, it appears that emission from 
the EUV and X-ray layers can only explain the lower and typical luminosity LVCs, but not the 
highest luminosity LVCs. 

It is instructive to estimate what physical conditions are required to produce [01] luminosities 
in the LVC as high as 10~^ — 10"^ Lq. Consider a layer on the disk surface of thickness i, 
temperature T, and extending to radius from which most of the [01] emanates. This top and 
bottom layer of the disk has hydrogen nucleus density n and vertical column A^. Because the [01] 
6300 A transition is AE/k = 23,000 K above the ground state, we require T> several thousand 
degrees K for significant emission. The emerging [01] 6300 A luminosity produced by the surface 
layers is then given 

L[oi] = 7rronen(0)7[o7]^, (35) 

where 7[o7] is the collisional excitation rate for electrons on atomic oxygen, n{0) is the density of 
atomic oxygen, and rie is the electron density. We account here for both the top and bottom of 
the disk, but recall that 1/2 of the luminosity (that directed to the midplane) is absorbed by the 
optically thick disk. Eq. (35) assumes to be less than the critical density Uecr- (HM89 give 
Hecr ~ 10^ cm~3 SO Tig < Tiger is generally satisfied). HM89 give 7[o/] = 8.5 x lO~^T2'^'^e~'^-^^'^'^ cm^ 
s~^, with T4 = T/W^ K. Oxygen rapidly charge exchanges with hydrogen and therefore at high 
temperatures (T >> 200 K) 0+/0 = H+/H. Therefore, and n{0) = x{H)no, where 

no — 3 X 10~^n is the gas phase density of oxygen in both O and O"*" and x{H) is the abundance 
of atomic hydrogen. It follows that 

L[oi] ^ 6.5 X W-^Xex{H)n5rl^N2oT^-^^e-''-^/^* Lq, (36) 
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where ri4 = To/IO^^ cm, 715 = cm ^, and A^20 = cm ^. This analytic exercise shows 

that to produce L^qi] ^ 10~^ — 10^ Lq in the LVC requires, for example, surface layers with n ~ 10^ 
cm-3, N ~ 10^0 cm-2, ~ 60 AU, T ~ 10"^ K, and Xe ~ 0.5. In the EUV layer Xe ~ 1, r4 ~ 1, 
and rf^ ~ 1 (recall most of emission arises from r ~ Vg). Therefore, the [01] luminosity from the 
EUV layer is approximately 

Lf^j( ~ 6.5 X 10-''x{H)rl^n5N2o Lq. (37) 

This equation shows the difficulty in producing the observed [01] emission from the EUV layer: the 
fraction of neutral gas x{H) is very low in the EUV layer, or equivalently, the fraction of atomic 
oxygen is very low. In addition, the emission mostly arises from r ~ ~ 10^^ cm, which is not 
large, and n^N2o rarely exceeds unity. 

On the other hand, in the X-ray layer the temperature is of order T/^ ~ 0.1 — 0.4 and x{H) ~ 1 
so that 

Lfo^^ ~ 6.5 X lO-^Xengrf^iVsoTO-^^e-^-sm l^. (38) 

The X-ray layer has Xe ^ 0.1, r^^ < 10, < 100 and A^20 ^ 10. Even inserting these upper limits, 
we find that the [01] luminosity is at most 10"^ Lq if T = 4000 K and 3 x 10"^ Lq if T = 2000 K. 
Therefore, the [01] luminosity is extremely sensitive to the gas temperature (and its variation in r 
and z) in the X-ray layer. The Meijerink et al (2008) model has high temperatures, ~ 3000 — 4000 
K, in the X-ray layer out to 10-20 AU in their case with relatively high X-ray luminosity of 2 x 10"^^ 
erg s~^. For this case, they therefore find L^j-^ ~ 5 x 10~^ Lq. Our model for the same case, 
however, gets temperatures in the X-ray layer in the range 1500 — 2500 K, and therefore we get an 
[01] luminosity of only 5 x 10"'^ Lq (see §4.4 below). However, keeping the X-ray luminosity the 
same but assuming a much softer power law spectrum [L{E) oc from 0.1 to 2 keV], we do find 
that the X-ray heated gas becomes hotter and the [01] luminosities approach 10""^ Lq. We discuss 
in §4.3 the reasons for the differences in temperature in our model relative to that of Meijerink et 
al. We conclude that it is unlikely that the X-ray layer can provide the highest [01] luminosities 
observed in some LVC sources, but that soft X-ray sources can produce the typical luminosities. 

The "transition layer" {xg ~ 0.5) between the fully ionized EUV layer and the partially ionized 
X-ray layer is also unlikely to produce either the highest observed [01] luminosities, or even the 
typical luminosities. The density in this layer is similar to the density at the base of the EUV layer, 
or n ~ 10^<^4(^ cm~'^ at r = = 7 AU (see Eq. 12). However, beyond Vg = 7 AU the density 
falls as r~^/^ so that most of the emission arises from r ~ 7 AU. In addition, the column of this 
transition layer is roughly the column for optical depth unity in EUV photons, or N ~ 10^® cm~^. 
Therefore, the small ro and low A^ conspire to produce only L^qi] 10~^ Lq. 

We plan to investigate the possibility that the source of the LVC arises from the shear layer 
produced when the protostellar wind strikes the surface of the disk obliquely and sets up an outward 
moving layer of shocked wind and entrained disk surface gas (a "shear" layer, see Matsuyama et 
al 2009), hereafter MJH09). MJH09 show that this layer can have columns N ~ 10^^ - 10^° cm'^ 
out to 100 AU. As noted earlier, the oblique wind shock {vs ^20 km s~^) is insufficient to ionize 
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hydrogen or helium to provide the electrons needed for [01] excitation. Therefore, we require the 
EUV and soft X-rays to ionize this layer. This shear layer is likely turbulent so that there may 
be rapid mixing of the bottom (cooler and more neutral layers) with the top shear layers, perhaps 
allowing Xe ^ 0.5 and x{H) ^ 0.5 in the entire layer (the most efficient for producing [01]) and 
maintaining a relatively isothermal layer. The heating would be a combination of shock/turbulent 
heating plus the heating due to photoionization by EUV and soft X-rays. The density n in the 
shear layer (MJH09) is approximately n ~ SOOOM-gr^g'^ cm~^. Note that we cannot allow M_8 to 
exceed unity or the EUV and soft X-rays will not be able to penetrate the base of the wind to heat 
and ionize the shear layer. Therefore, although the shear layer may provide sufficient column N , 
electron fraction x^-, and temperature T, it appears unlikely to produce sufficient nr^ to give the 
observed [01] luminosities in the more luminous sources. Further work is needed to confirm this 
rough argument. 

We conclude that the origin of the very luminous LVC [01] emission is not from the EUV layer, 
the X-ray layer, or the transition layer. The typical LVC [01] emission, however, may be produced 
by soft X-rays. The LVCs are also unlikely to originate from the shear layer set up by the impact 
of the protostellar wind. Perhaps a model that invokes ambipolar diffusion as a heating source, 
such as those that Safier (1993a, b) has constructed for the HVC, might be applicable for the most 
luminous LVCs. However, the EUV layer is capable of producing the lowest luminosity sources, 
and the X-ray layer may produce the typical luminosity, so we proceed with detailed numerical 
studies of the [01] luminosity from these layers in §4 below. 



4. NUMERICAL MODEL AND RESULTS 

4.1. The EUV Surface Layer and the Underlying X-ray Layer 

Gorti & Hollenbach (2008) describe the numerical code which we use to calculate, self-consistently, 
the gas temperature, gas density, and chemical structure of the predominantly neutral gas in the 
disk. To summarize, the code includes ~ 600 reactions among 84 species, gas heating by a number 
of mechanisms including FUV grain/PAH photoelectric heating and the heating caused by X-ray 
ionization ionization of the gas, and cooling not only from collisional excitation of the species fol- 
lowed by radiative decay, but also from gas- grain collisions when the dust is colder than the gas. 
In some instances, for example deep in the disk below the surface layers, the dust is warmer than 
the gas in which case gas-grain collisions can be an important heating source for the gas. The 
vertical structure of the disk is calculated self-consistently by using the computed gas temperature 
and density to calculate the thermal pressure, and then balancing thermal pressure gradients with 
the vertical (downward) gravitational force from the central star. 

For this paper we consider the "chemistry" of the fully ionized (HII) surface region, photoion- 
ized by the EUV radiation field from the star. By "chemistry" we mean the computation of the 
different ionic states of a given element by balancing photoionization with electronic recombina- 
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tion and charge exchange reactions. Photoionization rates are computed using cross sections from 
Verner et al. (1996). Recombination rates are taken from Aldrovandi & Pequignot (1973), Shuh &: 
van Steenberg (1982) and Arnaud & Rothenflug (1985). Charge exchange rates are from Kingdon 
&: Ferland (1996). We assume a gas temperature of 10^ K and do not perform thermal balance 
calculations in the HII region. In short, our code for the surface of the disk is an HII region code 
where we assume an EUV spectrum from the central star and then compute the abundances of, for 
example, Ne"*", Ne"'"'' and Ne"*"*"^ at each point r, z in the surface HII region. At each point r, z we 
compute both the direct EUV flux from the star and the diffuse EUV field caused by recombinations 
to the ground state of atomic hydrogen in the surface HII region. We use the method described by 
HoUenbach et al (1994) and utilized by Font et al (2004) to do both these computations. The code 
finds the electron density ne{r,zjF) at the base of the surface HII region (in other words, at the 
ionization front separating the ionized surface from the predominantly neutral disk below). The 
thermal pressure at the IF is then Pjp ~ 2n(,{r, ziF)kTij, where T// = 10^ K and the factor of 2 
includes the pressure from protons. He"*" and electrons. This pressure then determines the height 
zjF where the thermal pressure in the predominantly neutral gas has dropped from its midplane 
value to PjF- The parameter zjf is the height of the base of the HII region: above zif the emission 
is mostly EUV-induced and we call this region the EUV surface layer; below zjf the emission is 
mostly X-ray-induced and we call this region the X-ray layerQ 

Implicit in our model is the assumption that the EUV luminosity and the X-ray luminosity 
is generated close to the stellar surface. This assumption then allows us to determine the column 
density of wind that the EUV and X-ray fluxes must traverse (see §2) as well as the angle of 
incidence of the EUV and X-ray flux on the flared disk surface. Since we assume the protostellar 
disk wind originates near r '-^ 10^^ cm, our estimate of the attenuation column density at the wind 
base is valid as long as the EUV and X-ray source of luminosity originates roughly within this 
distance from the stellar surface. Models of X-ray flares indicate that the X-rays probably arise 
from flares whose size ranges from 0.1 to 10 times the stellar radius (e.g., Favata et al 2001, Grosso 
et al 2004, Favata et al 2005, Franciosini et al 2007, Stelzer et al 2007, Getman et al 2008a,b). 
Therefore, our estimate of the attenuation column is likely valid. Similarly, the line emission that 
we model usually arises from r -^1-10 AU in the disk, and so the placement of the EUV or X-ray 
source within 10^^ cm of the star does not affect our results. However, in DG Tau, a soft X-ray 
source has been imaged and seen to arise about 20 AU from the star, probably from shocks in a jet 
(Giidel et al 2008, Schneider & Schmitt 2008). Such a geometry would certainly lower the column 
of attenuating wind, because of the spherical divergence of the wind. In addition, the X-ray flux 
would strike the disk from above, nearly normal to the surface. This latter effect, however, will 
likely not significantly affect the luminosities in the lines, since as we have shown in §3, the line 
luminosities are really an emission measure effect, and mainly depend on the fraction of energetic 
photons that the disk absorbs. If the source is 20 AU from the star, roughly half of the energetic 



^We note that FUV may contribute to the heating in the X-ray layer. However, X-rays produce Ne^ and the 
electrons necessary for efficient excitation of [Nell] and [01]. Therefore, it is proper to call this the "X-ray layer". 
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photons are absorbed. We show below that in the case of a flared disk with a central source 
of energetic photons, nearly 0.7 of the photons are absorbed. Therefore the fraction of photons 
absorbed is nearly the same. 

Also implicit in our model is that the X-ray luminosity is the mean value of the time variable 
X-ray luminosity. Getman et al (2008a) show that typical decay times of flares is of order hours to 
days. In the EUV and X-ray layers where the modeled lines originate, recombination and cooling 
timescales are of order 1-10 years. Thus, the gas generally does not have time to respond to the 
flares, but settles to a state given by the mean value of the X-ray luminosity. 

4.2. The EUV Layer Results 

Figures 1 and 2 show the results of models as we vary the EUV luminosity. We assume two 
forms for the EUV spectrum. The first form (Figure 1) is a relatively hard spectrum; we assume a 
power law spectrum vL^ = constant from 13.6 eV to the X-ray regime (~ 0.1 keV). This spectrum 
is motivated by the fact that vLu in the EUV band is observed to be similar to vL^, in the X-ray 
band, and each band typically has vLy ~ 10~^ L*, where L* is the stellar bolometic luminosity. On 
the other hand, the EUV spectrum is very uncertain. The Ribas et al (2005) observations of older, 
but very nearby, solar mass stars show EUV spectra which can drop rapidly from 13.6 eV to 40 eV, 
even though the overall trend from the EUV to the X-ray tends to roughly an Ly oc z/~^ spectrum. 
To simulate a softer spectrum than the first form, we take a blackbody spectrum with an effective 
temperature of 30,000 K (Figure 2). We are further motivated to adopt a softer spectrum by the 
observations in one source of the ratio [Nelll] 15 //m/[NeII] 12.8 yum ^0.06 (Lahuis et al. 2007), 
and because [Nell] has been detected in more than 25 sources and none of them show [Nelll]. Our 
first form of EUV spectrum produces a ratio > 1! This either indicates very little production of 
[Nell] by the EUV layer, or that the EUV spectrum is much softer. We therefore have chosen a 
blackbody EUV spectrum that provides ratios in accord with measured values or upper limits on 
the ratio. 

Figures 1 and 2 show the nearly linear rise in I/jArg//] or L[Arem] with Lew as predicted in Eqs. 
(19) and (20). The absolute values are also in good agreement with these equations. At very high 
Lew this linear relationship breaks down for [Nelll], and Li\]eiii begin to saturate because the 
electron densities in the dominant emitting regions begin to exceed ni,cj.,[Neiii]- our power law 
spectrum in Figure 1 we find that by fitting our analytic results to the model that f{Ne~^^) ~ 0.75 
and f{Ne^) ^ 0.25; i.e., 75% of the emitting neon is in Ne"*""*" and only 25% in Ne"^. However, 
in Figure 2 we see that a softer EUV spectrum (blackbody with T(,ff = 30,000 K) will reverse 
the situation so that [Nell] dominates. Another mechanism, not treated here, that would quench 
[Nelll] and raise [Nell] would be turbulent mixing of neutral gas into the HII region. The charge 
exchange reaction Ne"''^ -|- H — > Ne"*" -|- H"*" is very rapid (e.g., Butler &; Dalgarno 1980), and even 
a neutral fraction x{H) ~ 10"^ would lead to [NeIII]/[NeII] < 1. 



- 23 - 



Figure 3 shows the results for a number of other fine structure transitions hsted in Table 1. 
One sees that [Aril] 7 //m is the strongest predicted line not yet observed. Again, the analytic 
estimates (Eq. 8 and Table 1) are very good. 

Figure 4 shows the radial origin of the EUV-induced emission in the [Nell], [Nelll], [Aril], 
and [SIII] lines. We use here as a standard case ^euv = 10^^ s~^ {Lew = 2 x lO^'' erg s~^) and 
the blackbody spectrum, and Lx = 2 x lO^'^ erg s~^ (with our standard X-ray spectrum, see §4.3) 
although the radial origin is quite insensitive to these parameters. We see that most of the emission 
arises from r ~ ~ 10 AU, as predicted in §3. We plot 47rr^ times the emergent flux from one side 
of the disk. This roughly gives the luminosity arising from both sides of the disk and from a region 
extending from 0.5r to 1.5r. For dominant ions like Ne"*", the luminosity scales as n'^H{r)r'^ tx r-"^/^ 
for r < Vg and rig < Uecr and as for Ue > Uecr- The luminosity scales as n'^r'^ for r > Vg and 
ne < n^cr, so that here it scales as (see Eq. 13), as seen in [Aril] and [Nell]. For non-dominant 
ions such as Ne"*"*", the scalings change because the fraction of neon in Ne"*"^ changes with radius. 

Hollenbach et al (1994) and Gorti & Hollenbach (2009) show that significant photoevaporation 
flows proceed in the 10^ K gas in the EUV layer at > 1 AU. Thus, although our models here are 
static, the emitting gas is actually rising (and rapidly turning radial, see Font et al 2004) at speeds 
of order the sound speed (or 10 km s~^) from the surface of the disk. As discussed in Hollenbach et 
al (1994), the electron density at the base is not much affected by this flow. Photoionization and 
recombination timescales are sufficiently short that steady state still applies to the computation 
of the ionization state of each element in the gas. Therefore, we expect our model results on the 
emitted luminosities in each fine structure line to be well approximated by the static model solution. 
However, the observed line profiles will be affected by this flow. For a disk viewed edge on, the 
lines will be broadened by the Keplerian rotation but also the radial outward flow. For a disk 
viewed face-on, the lines will be broadened mostly by the radial outward fiow, and since the far 
side of the disk is obscured, one would expect a blue shift. Alexander (2008) has recently modeled 
[Nell] line profiles from photoevaporating disks. He predicts broad (30 — 40km s~^), double-peaked 
profiles from edge-on disks due to rotation, and a narrower (~ 10 km s~^) profile with a significant 
blue shift (5 — 10km s~^) from face-on disks. He argues that the observed line widths in TW Hya 
(Herczeg et al. 2007) are consistent with a photoevaporative wind (see also Pascucci &: Sterzik 
2009). Resolved [Nell] observations can thus provide a test of EUV photoevaporation models. 

Figure 5 shows the vertical origin of the EUV-induced and X-ray-induced emission at r = 10 
AU for the standard case. We have plotted gas temperature T, the dust temperature T^ust, and 
the hydrogen nucleus density n as a function of the hydrogen nucleus column N measured from 
z = r (the putative "surface" of the disk) downwards. On the top of the figure we give the values 
of z that correspond to those of N. The completely ionized EUV layer extends to 3 x 10^*^ 
cm~^, and has electron densities ng ~ n ~ 3 x 10^ cm~^ (see Eq. 12). The X-ray heated (T ~ 1000 
K) and ionized layer extends from ~ 3 x 10^^ cm"^ to A ~ 3 x 10^'^ cm~^ with hydrogen atom 
densities ~ 3 x 10^ cm~^. FUV photons also contribute to the heating of this layer. 



4.3. The X-ray Layer 



We model the X-ray layer for both a "soft" X-ray spectrum and a "hard" X-ray spectrum. 
Our standard ("hard") X-ray spectrum (Gorti & Hollenbach 2004,2008) is based on observed X-ray 
spectra from young stars, with an attempt to correct for extinction at the softer energies (Feigelson 
& Montmerle 1999). Our fit to this spectrum is that of a power law (x v from 0.1 keV < hu < 2 
keV, fitting to another power law oc for hv > 2 keV. This spectrum is similar to that 
adopted by Glassgold et al (2007) and Meijerink et al (2008). We also model disks illuminated by 
a softer X-ray spectrum: Ly oc for 0.1 keV < hu < 2 keV, and Lj, oc z^"^'^^ for hv > 2 keVcl 
Ercolano et al 2009 recently provide evidence that such a spectrum might be expected from young, 
low mass stars. We note that our soft X-ray spectrum has equal energy flux in equal logarithmic 
intervals of photon energy between 0.1 keV and 2 keV; that is, there is as much energy flux between 
0.1 and 0.2 keV as there is between 1 and 2 keV. We do not consider here a harder spectrum than 
our "hard" case, although recently there have been observations of "superhot" flares (Getman et al 
2008a) that indicate significant emission in the 3-8 keV region of the spectrum. We do extend our 
"soft" and "hard" spectra to 10 keV, but there is insignificant energy flux beyond a few keV. If there 
were, then for the same X-ray luminosity as our two cases, we would obtain less emission in the 
lines we model in this paper. The higher energy photons penetrate more column of gas, depositing 
less energy per unit volume, and therefore lead to cooler gas than in our current X-ray layer. In 
addition, because of the increased penetration, the heat is deposited in molecular regions, where 
the cooling is enhanced by the molecular transitions. Therefore the emitting gas is substantially 
cooler and most of the X-ray heating energy presumably emerges in molecular rotation lines of, 
for example, CO, OH, and II2O or possibly, if grains are abundant, as IR continuum emission 
from grains heated by collisions with the warmer gas. However, we emphasize that the heating 
and cooling timescales are long, of order 1-10 years, so that superhot flares that are much more 
short-lived than this timescale will not produce a significant effect. 

Glassgold et al (2007) first modeled and discussed the physics of the X-ray layer, and our 
results are in basic accord with theirs, except, as discussed below, we obtain somewhat cooler 
temperatures in the X-ray layer. Figure 6 shows the vertical structure of the same case plotted in 
Figure 5 at the same radius, r = 10 AU, only the N range is shrunk to emphasize the X-ray layer. 
In Figure 6, we plot the electron abundance Xe relative to H nuclei, and the fraction of neon in Ne"*", 
fNe[Ne^). The [Nell] 12.8 ^Lm line luminosity is proportional to XefNe{^e~^), and very sensitive 
to T in the X-ray layer (see Figure 5 and Eq. (28), where the T dependence comes in the collisional 
rate coefficient, which is proportional to e~^^^^^^). Below the EUV layer to a depth A/" ~ 3 x 10^*^ 
cm~^. X-rays maintain a relatively high fraction of Ne"^, fNe{Ne~^) ^ 10~^, and X-ray ionization of 
H and He as well as FUV ionization of C maintain a relatively high ionization fraction, > 2 x 10^^ 
(see Figure 6). Note that the column attenuating the X-rays in this layer is the column through 



^ We note that our "soft" X-ray spectrum has the same power law L,^ <yiv ^ form as our "hard" EUV spectrum! 
In addition, we note that even our "hard" X-ray spectrum includes a contribution from 0.1 keV X-rays. 
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the disk to the star, which is typically ~ lOA'^, where N is the vertical column to the disk surface 
In the X-ray layer, 10^^ cm~^ < N < 10^^ cm~^, the gas is heated by a combination of FUV 
grain photoelectric heating and X-ray photoionization heating. It is cooled mainly by [01] 63 ^m, 
[01] 6300 A, [Nell] 12.8 fim, [Aril] 7 /im, and gas-grain collisions (see Gorti & Hollenbach 2008). 
The resultant temperature is ^ 1000 — 2000 K, dropping with increasing r and as dilution and 
attenuation of the X-rays and FUV lower the heating rates. Because T, fi^(,{Ne^), and Xg drop 
with increasing (see Fig. 6) and r, most of the [Nell] 12.8 /im emission arises from r < 20 AU 
and N ~ 10^° - lO^^ cm'^ where T ~ 1000 K, /(A^e+) ~ 10-\ and Xe ~ 10"^ 

In addition to showing the radial origin of the [Nell] emission in the EUV layer, Figure 4 also 
shows the [Nell] "luminosity" plotted as a function of r for the X-ray layer (hard spectrum X-rays) . 
Note that there is greater contribution from inner (~ 1 AU) regions of the disk compared to the 
EUV layer. In addition, there is more luminosity emerging from the X-ray layer than the EUV 
layer. Figure 4 shows that the X-ray-induced emission arises mostly from r ~ 1 — 10 AU, also in 
agreement with Glassgold et al. Beyond this radius, the X-ray and FUV heating is insufficient to 
maintain significant quantities of T > 1000 K gas. 

Figure 7 plots the [Nell] 12.8 /im and [Nelll] 15 /xm emission from the X-ray layer. In agreement 
with Glassgold et al (2007) and Meijerink et al (2008) we find that [Nelll] / [Nell] ^0.1, mainly 
caused by the rapid charge exchange of atomic H with Ne"^"*" . We also plot outlines of the observed 
54 sources tabulated by Giidel et al (2009). The vertical dotted lines shade the region where the 
sources have known outflows and jets. The horizontal solid lines shade the regions with no outflows 
or jets detected. The Lx tabulated by Giidel et al is a two component fit with an attempt to 
correct for absorption of the softer X-rays by material on the line of sight from star to observer. 
However, many of the observations do not extend to hv < 0.3 keV and extinction is severe at 
the lower energies, so that a luminous soft X-ray source that is weak at 0.5-1 keV could exist 
undetected. The effect of such a "soft" X-ray component would be to move the data points to the 
right on Figure 7 and comparison should be made to our "soft" X-ray spectrum results (dashed 
line). We also find that Lj^^ii and LncIII scale with Lx-, as predicted in §3 and also as found by 
Meijerink et al (2008). Comparison with Figure 2 shows that if the X-ray luminosity is about the 
same as the EUV luminosity from the central star, and if the EUV spectrum is soft enough that 
[Nell] dominates [Nelll] in the EUV layer, the [Nell] luminosity is roughly 2 times stronger from 
the X-ray layer as from the EUV layer, as we estimated analytically. The main conclusion from 
comparing the data to the model results is that although the X-ray layer may explain the origin of 
the [Nell] emission in some (perhaps most if a strong soft X-ray excess is common) sources, there 
are a significant number of sources, especially those with observed outflows and jets, where the 
X-ray luminosity seems insufficient to explain the [Nell] luminosity. In §4.4 below we compare the 
observational data with our analytic results on the [Nell] luminosity expected from internal shocks 
in the jets or winds, and find that this is a plausible origin for these sources. 

As noted above, our results on the IR fine structure emission from the X-ray layer do not differ 
appreciably from previous results (Glassgold et al 2007, Meijerink et al 2008). Overall, we tend to 
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find somewhat (factor of ^ 2) lower IR line luminosities. This agreement is a bit fortuitous, arising 
from a cancellation of several effects and the insensitivity of the fine structure lines with variations 
in T if T > 1000 K. Our models self-consistently calculate the vertical density structure of the gas 
by using the computed gas temperatures (which differ from the dust temperature) to calculate the 
gas density structure rather than assuming the gas density structure is fixed by the calculation of 
vertical pressure balance when one assumes the gas temperature to equal the dust temperature, 
as done in the previous work. Our self-consistent model produces significantly different results at 
columns A'^<10^^ cm~^, where the gas temperature rises above the dust temperature (Gorti & 
Hollenbach 2008). The net effect is that our gas disk is more fiared, intercepting a larger fraction 
of the X-ray luminosity. This tends to raise the emission from our models. In addition, we include 
FUV grain photoelectric heating which also raises the emission. However, counteracting these 
effects is the inclusion of more gas coolants in our model, especially [Nell] 12.8 fj,m and [Aril] 7 
/im. In addition, our treatment of the gas heating by X-rays follows Maloney et al (1996), which is 
somewhat different than the approach used by Glassgold et al and Meijerink et al, and our X-ray 
heating rates are lower than these authors by a factor 3-10. We believe this may arise because we 
include the loss of " heat" due to escape of Lyman alpha and other photons created by recombining 
hydrogen or to the absorption of these photons by dust. Overall, our X-ray layer tends to be a 
factor of about 2 cooler than the previous models (roughly 1000-2000 K versus 2000-4000 K in 
the previous models), thereby lowering the fine structure emission from this layer. This lower 
temperature, has a relatively small effect on the fine structure lines, because their upper states lie 
only AE/k ~ 1000 K above the ground state. However, it has an enormous effect on our predictions 
of the [01] 6300 A emission, whose upper state lies AE/k ~ 23,000 K above the ground state, as 
we will discuss below in S4.5. 



4.4. [Nell] Emission from Internal Shocks in the Jets and Winds 

Figure 8 plots the [Nell] luminosity versus the mass accretion rates assembled by Giidel et al 
(2009). As in Figure 7, the vertical dotted lines shade the region which includes sources with known 
jets or outflows, whereas the solid horizontal lines denote sources with no detected jets/outfiows. We 
plot here our predicted [Nell] luminosities from internal shocks in the winds/jets, using our analytic 
expression (Eq. 33). The solid line represents the expected [Nell] luminosity when = O.lMacc, 
the entire wind or jet passes through a shock or /^/j = 1, the shock velocity is in excess of about 100 
km s~^, and the preshock density is less than 10'^ cm~^. The upper dashed line makes the same 
assumptions except that = Mace and the lower dashed line assumes = O.OlMacc- Note that 
the [Nell] luminosity is proportional to the product of fsh and M^, so that, for example, the lower 
dashed line also corresponds to fsh = 0.1 and = O.lMacc- The main conclusion is that internal 
wind or jet shocks very likely explain the origin of the [Nell] from the outflow and jet sources. In 
fact, the flgure might suggest that these shocks could explain the [Nell] observed in nearly all of 
the sources, if it were not for the fact that in some cases (e.g., Herczeg et al 2007, Najita et al 
2009, Pascucci et al 2009) where the lines have been spectrally resolved, they are narrower than a 
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shock origin would predict. We do note that in some of these cases, the integrated flux seen with 
the high spectral and spatial resolution ground based instruments is significantly less than the flux 
seen by the low resolution Spitzer Space Telescope. Najita et al speculate that perhaps there are 
two components comprising the total flux, a strong but broad and extended shock component and 
a weaker, but narrow and spatially unresolved disk component arising from the X-ray layer. On 
the other hand, it is quite possible that X-rays or EUV produce most of the [Nell] luminosity in 
the sources with no observed winds or jets. Note that these sources in Figure 8 are distributed in a 
nearly horizonal line with no apparent dependence on Mace over two orders of magnitude increase 
in this parameter. 

4.5. [OI] 6300 A Emission from the EUV and X-ray Layer 

Figure 9 plots the [01] 6300 A luminosity from the EUV layer versus ^euv for both our harder 
LEUvi^) spectrum and our softer Leuv{i^) blackbody spectrum. A harder spectrum gives 

more [01] luminosity in the EUV layer because, although the gas is almost entirely ionized, there 
is a greater fraction of neutral H and O in the gas due to the smaller photoionization cross section 
of these atoms with higher photon frequency. However, even the harder EUV spectrum results in 
[01] 6300 A luminosities < 10"^ Lq, which can only explain the weakest LVC sources. Recall that 
L[oi] ranges from 10"^ — 10"^ Lq in LVCs. 

Figure 10 plots the [01] 6300 A luminosity in the X-ray layer versus Lx for both our harder 
and our softer X-ray spectrum. The [01] luminosity does increase with Lx, due to the higher 
temperatures and higher ionization fractions in the X-ray layer. Recall that the gas is primarily 
neutral, so the higher ionization fraction increases the luminosity by increasing the density of the 
electrons which excite the [01]. However, with our standard (harder) X-ray spectrum, which is quite 
similar to that adopted by Meijerink et al (2008), we obtain [01] luminosities that are a factor of 
nearly 100 times lower than those of Meijerink et al (2008) from the X-ray layer. This is primarily 
because of the extreme sensitivity of the [01] luminosity to the temperature of the X-ray layer (see 
Eq. 38). As noted above, our temperatures are roughly a factor of 2 lower than those of Meijerink 
et al. 

The main point, however, is that in our standard models neither the EUV layer or the X- 
ray layer can produce [01] luminosities as high as 10~^ — 10~^ L0 as observed in many LVC 
sources. This confirms the analytic estimates made in Section 3.5. However, as discussed above, 
the [01] luminosity is very sensitive to the temperature of the X-ray layer. One way to increase the 
temperature is to assume a softer X-ray spectrum. Softer X-rays have much higher absorption cross 
sections, and therefore deposit much more heat per unit volume in the upper layers. In addition, 
they create higher electron abundances, and these lead to increased efficiency in converting the 
absorbed X-ray energy into heat. Therefore, we also plot in Figure 10 the results for cases with 
similar X-ray luminosities, but with our "soft" X-ray spectrum where L^, oc from 0.1 keV to 
2 keV. This spectrum has many more 0.1-0.3 keV X-rays than our standard case, and we find 
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that we do indeed get higher temperatures and electron abundances in the upper parts of the 
X-ray layer, and consequently much higher [01] luminosities, -^^[o/] can be as high as ~ 10~^ 
Lq if Lx ~ 10^^ erg s~^ (~ 2 x 10~^ Lq), a likely upper limit to the soft X-ray luminosity. 
Therefore, soft X-rays may be able to explain "typical", L^q^ ~ 10~^ L©, LVC sources, but not 
the most luminous sources. We note that if the soft X-rays were causing photoevaporation, then 
high spectral resolution observations of the [01] line might diagnose the flow parameters. 

5. DISCUSSION 

This paper has focused on fine structure lines from ions which required hv > 13.6 eV photons 
to photoionize them. Most of these lines fall in the 5 /im < A < 40 /im wavelength region and are 
therefore partially accessible through atmospheric windows from the ground, and entirely accessible 
from space-based observatories such as the Spitzer Space Telescope. Ground based observatories 
have the advantage of larger diameter telescopes and therefore greater spatial resolution as well as 
larger and heavier instruments capable of higher spectral resolution. The TEXES instrument (Lacy 
et al. 2002) achieves a spectral resolution of ~ 3 km s~^ and a spatial resolution of ~ 0.5(A/10;um) 
arc seconds on a 10 meter-class ground-based telescope such as Gemini. Its sensitivity to line flux 
(5 sigma in one hour) translates to L ~ 3 x 10^^ Lq at 100 pc. The Michelle instrument (Glasse 
et al. 1997) achieves spectral resolution of ~ 15 km s~^ and is capable of detecting lines with 
luminosities L ~ 3 x 10~^ at 100 pc, if mounted on a 10 meter-class ground-based telescope. 
The VISIR instrument on an 8 meter class telescope has a sensitivity of ~ 3 x 10~^ Lq at 100 
pc, and a spectral resolution of about 12 km s~^ (Lagage et al 2004). Spitzer had relatively poor 
spatial (~ 12 arcsec resolution ) and spectral 500 km s~^) resolution but could achieve 5 sigma 
in one hour sensitivity that translated to L '--^ 10~^ Lq at 100 pc. 

A number of groups have now observed nearby star-disk systems and measured fluxes from 
especially the [Nell] 12.8 fiui line, with a few detections of H I reco mbination lines and good upper 
limits for [Nelll] 15 /im lines. Many of the observations were done using the IRS spectrometer on 
Spitzer (Espaillat et al 2007, Lahuis et al 2007, Pascucci et al 2007, Ratzka et al 2007). 

However, we first discuss recent ground-based observations with high resolution spectroscopy 
that help constrain the origin of the Nell emission by interpretation of the observed linewidths and 
spatial extents, as well as by the observed fluxes (e.g., Herczeg et al. 2007, van Boekel et el. 2009, 
Najita et al. 2009, Pascucci & Sterzik 2009). The first such resolved source to be observed and 
also one of the brightest is TW Hya (Herczeg et al. 2007). Herczeg et al interpret the observed line 
width (~ 21 km s~^) from this nearly face-on disk as possibly indicating the emission arises from 
the inner regions (~ 0.1 AU) of the disk. In our models it is very difficult to produce the observed 
[Nell] luminosity from X-rays or EUV at 0.1 AU. However, as they note, it might also originate 
from the EUV or X-ray layers at r ~ 10 AU, if turbulence can produce the observed linewidths. 
Alternatively, the linewidth may arise from the fact that the gas is not merely in Keplerian rotation, 
but also is photoevaporating at 10 km s~^ with respect to the disk surface. This produces a 
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blue shift of the [Nell] with respect to the stellar velocity (e.g., Alexander 2008). Using the VISIR 
spectrograph on the VLT telescope Melipan, Pascucci &; Sterzik (2009) recently showed that nearly 
all the line flux is blueshifted, with a peak at -6.3 km s~^ and a FWHM of 14.2 km s~^. They point 
out that these observations are in near perfect agreement with the prediction of Alexander (2008) 
for [Nell] produced in EUV-induced photoevaporating flow, and inconsistent with a static disk 
atmosphere. Alternatively, a soft X-ray spectrum might produce a very similar photoevaporating 
profile, since soft X-rays heat the disk surface at 1-10 AU to almost the same temperatures as 
the EUV layer. We note that TWHya is known to have a strong soft X-ray excess (Kastner et 
al 2002). The measured low accretion rate, 5 x 10^^*^ Mq yr~^ (Muzerolle et al. 2000), and the 
absence of any known outflow support an EUV and/or X-ray heated disk origin for the [Nell] 
emission. The observed [Nell] luminosity as measured by Pascucci &; Sterzik is ~ 4 x 10~^ Lq. 
From our models, we predict no contribution from shocks which is consistent with the low observed 
linewidths. We would expect EUV and X-rays to irradiate the disk given the low accretion rate and 
the measured flux to be a sum of the contributions from the ionized and neutral layers of the disk. 
We calculate the contribution from the neutral layer to be ~ 3 x 10~^ Lq (using Lx = 2 x 10'^'^ erg 
s~^; Kastner et al. 2002). If the remaining ~ 10~^ Lq is from the ionized layer, we then estimate 
that ^EUV = 3 x 10^*^ s~^ for TW Hya. However, given the accuracy of our models, we cannot 
rule out that most of the [Nell] emission is from the EUV layer {^euv = 1-2 x 10"^^ s^"*^) . The 
excellent agreement of the [Nell] line profile with the EUV model suggests that EUV may dominate 
in this source, but our model results using the observed Xray luminosity suggests that a substantial 
amount of the [Nell] may arise in the X-ray layer. Modeling of X-ray induced flows and further 
observations are needed to clarify this discrepancy, possibly of the [Aril] 7 /im line which might 
discriminate between [Nell] emission from EUV or X-ray layers (see discussion below at end of this 
section). 

Herczeg et al. fail to detect [Nell] emission for the sources BP Tau and DP Tau. These non- 
detections are also compatible with the [Nell] emission models described in this paper. BP Tau is 
an actively accreting star ( ~ 2 x 10~^ Mq yr~^; Muzerolle et al. 2000) with presumably no EUV 
penetration of the disk wind and an X-ray luminosity {Lx ^ 7 x 10^^ erg s~^) that would produce 
lower [Nell] emission than the upper limit from observations. DP Tau has a low accretion rate, but 
very poor upper limits to the [Nell] flux to provide any reasonable estimates of ^euv- 

Van Boekel et al. (2009) report that the [Nell] emission from the T Tau triplet, which is 
resolved spatially and spectrally, has large linewidths ~ 100 km s~^ and is associated with a known 
outflow. Even though T Tau N is a very strong X-ray source with Lx ~ 2 x 10'^^ erg s~^ (Giidel 
et al. 2007), the expected [Nell] from the disk is still a factor of ~ 10 lower than what is observed. 
In addition, van Boekel et al spatially resolve the [Nell] emission and determine that it arises from 
T Tau S. We do not expect EUV and soft X-rays to penetrate the disk wind for this young source. 
On the other hand, our models of shock emission are consistent with the van Boekel et al data, 
as discussed in §3.4. van Boekel et al also conclude that shock emission is the likely origin of the 
[Nell] emission. 
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Najita et al. (2009) have observed two young disks around AA Tau and GM Aur using TEXES 
on Gemini N, and spectrally resolved the [Nell] line in both the sources. The FWHM linewidths 
are 70 and 14 km s~^, respectively, and the authors interpret the emission as arising from the X-ray 
layer in Keplerian disks. They also note that the flux in the line is less than that measured by the 
much larger beam of Spitzer. A spatially extended and broad (FWHM) additional component, such 
as a protostellar wind shock, could account for the difference. GM Aur is a transition disk object, 
that is still actively accreting at ~ 10~^MQyr~^, indicating the presence of gas in the dust depleted 
inner disk. The disk accretion rate is at a marginal epoch where the EUV may make it through to 
irradiate the disk or may be absorbed by the disk wind. From the known X-ray luminosity of the 
star {Lx ~ lO^^erg s'^; Strom et al. 1990), we estimate an X-ray produced [Nell] line luminosity 
of 2 X 10"^ Lq, while the observed value is ~ 7 X 10-6 Lq (Najita et al. 2009). The rest may arise 
from shocks, although no known outflows exist. Alternately it may come from either an unobserved 
EUV or soft X-ray component that has just begun to penetrate the disk wind and heat and ionize 
the surface layers. Note that if EUV dominates , ^euv ~ 2 x 10^^ s"^. The classical T Tauri star, 
AA Tau has a low accretion rate for an object of its class, estimated at 3 x 10~^M.q yr"^ (Gullbring 
et al. 1998), and we expect irradiation of the disk by EUV and X-ray photons due to the expected 
low wind column density. The observed line luminosity is ~ 4 x 10"^ Lq. AA Tau is highly X-ray 
variable with Lx ~ 3 X 10^° - 2 X 10^^ erg s"^ (Schmitt & Robrade 2007) which can result in Nell 
luminosities arising from the X-ray heated neutral layer, ranging from lO-^L© to IO'^Lq, and the 
observed value lies within this range. While it is likely that the X-ray layer explains the origin of 
the [Nell], the observed [Nell] flux places an upper limit of ^euv ^ 2 x 10^^ for AA Tau. 

Pascucci &: Sterzik (2009) detect [Nell] in all the three transition disks that they observe (TW 
Hya, CS Cha, T Cha), but only from one of the three classical disks (Sz 73). They claim that the 
resolved linewidths of all the transition disks are consistent with a photoevaporative flow driven 
by stellar EUV photons and estimate ^euv ~ lo^^"^^ s"^. These numbers should be considered 
as upper limits to ^euv since there may be some contribution to the [Nell] flux from the X-ray 
layer. Pascucci &: Sterzik also observe blueshifted [Nell] emission in CS Cha and T Cha, consistent 
with EUV photoevaporation. In CS Cha the inferred hole size is 45 AU. We note that if the 
inner disk is completely clear of gas, such a large hole is only consistent with EUV- induced [Nell] 
emission since the X-ray flux at this radius is too low to heat the gas to temperatures > 1000 K 
required to excite the [Nell]. Pascucci & Sterzik point out that they only detect such evidence of 
EUV photoevaporation in sources with very low accretion rates, consistent with our model here 
that EUV does not penetrate the wind base until the accretion rates are low. We also note that 
the expected X-ray heated [Nell] emissons for their non-detections are consistent with their upper 
limits. 

We next discuss the totality of [Nell] observations, which is dominated by unresolved Spitzer 
sources. There is clearly a considerable amount of scatter when one tries to see if the [Nell] 
luminosities Lj^eii correlate with either the X-ray luminosity Lx or with Mace- Espaillat et al. 
(2007) conclude that the [Nell] has a nearly linear correlation with the mass accretion rate; they 
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find a 10 times increase in [Nell] luminosity with about a 10-fold increase in accretion rateO 
However, Espaillat et al had a data set of only 7 sources, whereas recently, Giidel et al (2009) have 
compiled a data set of more than 50. Giidel et al find little L]\ren dependence with mass accretion 
rates at low Mace ^ 3 x 10~^ M0 yr~^, but a roughly hnear trend in the high Mace ^ 3 x 10~^ M0 
yr~^ sources which show evidence for jets and outflows (see Figures 7 and 8). The latter suggests 
that protostellar wind shocks may be responsible for the [Nell] from the outflow sources. It seems 
unlikely that the [Nell] in these sources is due to soft X-rays or EUV, since the wind mass loss rates 
are sufficiently high to likely block these photons from ever striking the disk surface at radii near 
rg. In addition, our analytic predictions of [Nell] luminosities from wind shocks seem to match 
the observations (Figure 8). Espaillat et al find little correlation of LmcII with Lx- Giidel et al 
formally find L^eii oc L^^^ but with a tremendous amount of scatter. We note that although many 
of the observed X-ray luminosities derive from observations of ~ 0.2 — 10 keV X-rays, the soft (0.1- 
0.3 keV) X-rays may suffer considerable extinction that is difficult to estimate, and considerable 
luminosity could be "hidden" in such a soft component. Some of the observed scatter may then be 
caused by [Nell] arising from EUV, soft X-ray, or shock heated and ionized gas. 

In summary, shocks may dominate at high Mace ^ 3 x 10~^ Mq yr~^, but there is observational 
evidence that EUV or X-rays must dominate at lower accretion rates. Because X-rays are more 
efficient in producing [Nell], naturally produce [Nell] stronger than [Nelll] as observed, and more 
easily penetrate the base of the protostellar wind, it seems likely that X-rays often dominate the 
EUV production of [Nell] in disks, although not by a large factor. A part of this evidence for 
a non-shock origin has been gathered by high spectral resolution observations of [Nell] made by 
ground-based telescopes, which show relatively small linewidths compared to the ^ 100 km s~^ 
linewidths expected for wind-shocked [Nell]. Although ~ 1 keV X-rays may play a role in the 
production of L]\jeii for sources with weak winds, there is clear evidence that EUV or soft X-rays 
may sometime dominate. If one wanted to identify a source where it is likely that either EUV 
or soft X-rays dominate the [Nell] production, one would choose sources with low accretion rates, 
Mace ^8 x 10~^ Mq yr~^, whose LncII fies well above the observed correlation of L^eii with the 

1 keV Lx. 

One of our principal results is that X-rays are more efficient in producing [Nell] emission than 
are EUV photons. If the central star has the same luminosity in X-rays as it does in EUV photons, 
the [Nell] luminosity from the X-ray layer will be about 2 times greater than the [Nell] from the 
EUV layer (assuming a soft EUV spectrum, which is the most efficient in producing [Nell]). This 
result was shown both analytically, in Section 3, and in our numerical results as seen in Figures 

2 and 7. Since the [Nell] luminosity scales linearly with the EUV luminosity and with the X-ray 
luminosity, this means that the EUV luminosity needs to be at least 2 times the X-ray luminosity 
for the EUV to dominate the production of [Nell]. Unfortunately, we have little idea of the EUV 



^We note that Pascucci et al (2007) found a tentative anticorrelation with accretion rate, but this was based on 
a very limited dataset which had a small range in values of line luminosity and accretion rates. 



- 32 - 



luminosity, since it is impossible to observe in young sources. Observations of older, nearby stars 
by Ribas et al (2005) suggest that the luminosity in the EUV band is usually similar to that of the 
X-ray band. However, these sources are not accreting, and it is possible for an accreting source to 
be very bright in EUV relative to 1 keV X-rays (but see Alexander 2004b and Glassgold et al 2009, 
who argue these photons are attenuated by the accretion columns near the star). In any event, 
these accreting stars need to have sufficiently low wind mass loss rates to allow these accretion 
shock-generated EUV photons to penetrate the wind base and strike the outer disk to create EUV- 
generated [Nell]. Alexander et al (2005) estimate EUV fluxes from stars with observed ultraviolet 
emission lines, and conclude that in some cases the EUV photon luminosities can be as high as 
10^^ s~^. This suggests that in some cases, the chromospheric emission may generate more EUV 
luminosity than X-ray luminosity in young stars. However, our own results place upper limits on 
the possible EUV photon luminosities: ^euv ^ 10^^ s~^. Overall, it appears that it is unlikely that 
the EUV fluxes on the disk surface are any stronger than the X-ray fluxes, and that it is likely that 
X-rays often slightly dominate EUV photons in the production of [NeH] when the wind mass loss 
rates are low so that internal wind shocks are weak. 

We have plotted the observed [NeH] and [NeHI] data on Figure 7, using the compilation of 
Giidel et al (2009) which uniformly treats all previously observed sources. The observed [NeHI]/[NeH] 
ratio of less than 0.06 in the source with measurements of both lines (Szl02, Lahuis et al. 2007) 
favors either an origin in the X-ray layer, a shock, or in a soft {T^ff < 40,000 K) EUV layer. A 
hard EUV layer such as our adopted power law oc is ruled out. Several of the sources have 
[NeH] luminosities readily explained as arising in the X-ray regions, as noted by Meijerink et al 
(2008). However, a number of the sources have larger [NeH] luminosities than can be explained 
by ^ 0.5 keV X-rays alone. In many such cases, such as the T Tau South source discussed by van 
Boekel et al (2009), shocks in the protostellar wind are the likely source. We note that since wind 
mass loss rates scale with accretion rates, shocks would provide the observed correlation (Espaillat 
et al 2007, Giidel et al 2009) between the mass accretion rate and i[Are//] (see Figure 8). 

If the wind mass loss rates are not sufficient to provide the observed [Nell] luminosity, or if 
ground-based observations reveal narrower lines than might be expected from the shocks, such as 
in TW Hya, [Nell] emission may be generated by a soft EUV or X-ray spectrum from the central 
star. Since the "hard" X-rays were insufficient to explain some of these sources, and since we 
have shown that X-rays are more efficient in producing the [Nell] line, the only way that EUV 
luminosity from the central star can explain these sources is for the EUV luminosity to be greater 
(> 2 times) than the observed X-ray luminosity, and, in addition, the EUV spectrum has to be 
"soft" (Tg// < 40,000 K). If EUV does dominate, we can see from Figures 1 and 7 that Lew ^ 10"^ 
L0 is often required. A luminosity of 10^^ Lq corresponds to ^euv ~ 10^^ s~^. The comparison 
of the [Nell] and [Nelll] data with Figures 1 and 2 gives hard upper limits on ^euv- Most sources 
have ^EUV ^ lO'^^ s~^. If EUV is the main excitation mechanism, the comparison actually measures 
^EUV and the [NeII]/[NeIII] ratio constrains the EUV spectrum. 

We have examined our model results for diagnostics that would reveal whether the [Nell] 
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emission arises from the EUV layer or from the X-ray layer. One possible diagnostic is the ratio 
of the [Nell] 12.8 fim line to the [Aril] 7 fim line, [Nell] /[Aril]. We have shown in Figures 2 and 
3, along with the analytic calculation Eq. (8) combined with Table 1, that in the EUV layer the 
[Nell] / [Aril] ratio is about unity for our "soft" EUV spectrum. This spectrum maintains most Ne 
and Ar in singly ionized form in the EUV layer. Although elemental Ar is 20 times less abundant in 
the HII gas than Ne, the rate coefficient for electronic excitation of the [Aril] line is about 10 times 
larger than that of [Nell], and the 7 /um line has almost twice the photon energy as the 12.8 /um line, 
making up for the abundance discrepancy. In the X-ray layer, most of the Ne and Ar is neutral, 
and the fractional abundance of Ne+ and Ar"*" depends, in addition to elemental abundances, 
on the X-ray photoionization cross sections of Ar and Ne, on the electron rate coefficients for 
collisional ionization of Ar and Ne by secondary electrons, and on the rate coefficients for electronic 
recombination of Ne^ and Ar"*". In addition, the [Aril] line lies AE/k ~ 2060 K above ground, 
whereas the [Nell] line lies only ~ 1100 K above ground. Since the X-ray heated gas is typically 
~ 1000 K, this means that the relative line strengths are sensitive to the temperature of the X-ray 
heated layer, with [Aril] gaining advantage in warmer gas relative to [Nell]. We find in our models 
that for our hard X-ray spectrum, which peaks at 2 keV and where the Ne and Ar are ionized mainly 
by direct X-ray photoionization, the X-ray layer produces [Nell]/ [Aril] ~ 2.5. Unfortunately, due 
to a coincidence of atomic parameters and the enhanced heating due to soft X-rays, for our soft 
X-ray spectrum the ratio is [Nell]/ [Aril] ~ 1, the same as in the EUV layer. Thus, this ratio may 
discriminate between [Nell] produced in the X-ray layer versus the EUV layer only when the X-ray 
spectrum is relatively "hard". Nevertheless, a large ratio would strongly point to an origin in the 
X-ray layer. 

We have also examined both analytically and numerically the expected [01] 6300 A luminosity 
from disks around young stars. The observed luminosities in this line range from 10"^ to 10-3 Lq 
in the low velocity component, which has been identified as arising from the disk. We have shown 
that the EUV, transition, and (hard) X-ray layers are not likely to produce [01] 6300 A luminosities 
greater than 10"^ Lq. Meijerink et al (2008) provided models utilizing a relatively "hard" X-ray 
spectrum (peaking around 1 keV) which achieved [01] luminosities as high as ^ 10~^ L©, but our 
models with a similar X-ray spectrum give [01] luminosities ~ 10"^ Lq. The [01] 6300 A line is 
extremely sensitive to the temperature in the X-ray layer, as we showed analytically in Section 
3.5. Our models give typical temperatures of 1000 - 2000 K, whereas the Meijerink et al models 
give 2000-4000 K. We discussed in §4.3 the improvements in our models which lead to lower gas 
temperatures in the X-ray layer. However, Ercolano et al (2009) appeal to observational constraints 
on the emission measure distribution as a function of temperature for the chromospheres of young 
star analogs to argue that there is a (largely unobserved) soft X-ray component that is much larger 
than that assumed in our standard X-ray spectrum and in Meijerink et al (2008). Ercolano et al 
find that the X-ray spectrum may be better approximated by a power law L^j oc z^~^ from 0.1 keV 
to 2 keV. We have also run cases with such a soft X-ray spectrum, and find that X-ray luminosities 
of ~ 1032 erg s ^ can then give rise to [01] luminosities of ~ 10 ^ Lq. 
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6. SUMMARY AND CONCLUSIONS 

Circumstellar disks around low mass stars evolve with time with a decreasing accretion rate 
onto the star and a decreasing wind mass loss rate from the inner disk. X-rays, EUV and FUV 
photons from young, low mass stars arise principally from either magnetic activity (an active 
chromosphere) or from the accretion shock arising as disk material falls onto the star, presumably 
in accretion columns along stellar magnetic field lines. In the latter case, the energetic photons 
must penetrate or obliquely avoid the accretion columns in order to illuminate the disk surface. 
In either case, they must penetrate the protostellar wind near the wind base. We treat here the 
penetration of the protostellar wind and find that FUV photons likely penetrate first, when the 
wind mass loss rate is ^ 4 x 10~^ Mq yr~^, the exact number depending on the very uncertain 
dust opacity in the wind base material. As the wind mass loss rate drops with time, ~ 1 keV 
X-rays penetrate next, when Myj ~ 4 x 10~^ Mq yr~^. Finally, soft (~ 0.1 keV) X-rays and EUV 
photons penetrate only when the wind can be fully ionized at the base, which occurs roughly at 

^ 8 X 10"^'' Mq yr~^. The corresponding mass accretion rates onto the star are about 10 times 
higher, with considerable scatter. Considering observed rates of mass accretion with time (e.g., 
Hartmann et al 1998), these criteria translate to FUV and 1 keV X-rays penetrating very quickly 
after mass infall onto the disk from the molecular core has ceased, whereas EUV and soft X-rays 
may require an additional 1-2 Myr (with a lot of scatter) before they illuminate the disk. 

The 1 keV Xrays and FUV photons penetrate the disk surface to vertical columns of A'^ ~ 10^^ 
cm~^, and heat this layer to temperatures of order 1000 K for r < 10 — 20 AU. The X-rays ionize 
hydrogen and atoms with ionization potentials (IP) > 13.6 eV in this predominantly neutral layer, 
providing both electrons and species such as Ne"*" and Ar+. Thermal collisions of the electrons with 
these species produce fine structure lines such as [Nell] 12.8 ^m. The high gas temperatures and 
elevated electron abundances also produce strong emission from the [01] 6300 A forbidden line in 
regions with T>2000 K. The FUV photodissociates molecules, ionizes species with IP < 13.6 eV, 
and contributes to the gas heating. 

The EUV photons incident upon the disk create a fully ionized (HII) layer with T ^ 10^ K, 
that lies above the X-ray layer on the disk surface. Here, EUV photoionizes species with IP ^ 13.6 
eV, and singly or doubly ionized species tend to be the dominant ionization stage. Trace amounts 
of atomic oxygen are present and a relatively small amount of [01] 6300 A luminosity emerges from 
this layer. Due to a combination of falling electron density, rising scale height, and increasing disk 
surface area with increasing r, most of the fine structure emission from the EUV layer arises from 
r ~ ~ 7(M*/1 Mq) AU. The EUV layer produces more hydrogen recombination line luminosity 
than the X-ray layer, but does not explain the observed high ratio of these lines to [Nell]. It is 
likely the hydrogen recombination lines are produced in dense plasma close to the star: in the 
chromosphere, the accretion shock, or in a wind shock very close to the star. 

Strong ( > 100 km s~^) shocks, such as can be produced in internal shocks in protostellar winds 
or jets, can also significantly ionize species with IP > 13.6 eV and heat the gas to T >> 1000 K, 
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sufficient to excite the fine structure lines, the hydrogen recombination fines, and optical forbidden 
lines such as [01] 6300 A. Such ionization and heating has been inferred by observation of optical 
lines emitted in knots in the jets and in Herbig-Haro objects. 

In this paper we have analytically modeled all three of these emitting regions, and have pre- 
sented results from detailed thermo/chemical numerical models of the EUV and X-ray layer. We 
have focussed on the emergent line luminosities of [Nell] 12.8 /im, [Nelll] 15.5 ^m, [Aril] 7 /im, 
[Arlll] 9 /im, [SIII] 19 fiui, [SIII] 33 /xm, and [01] 6300 A. However, we also discussed infrared hy- 
drogen recombination lines (6-5 and 7-6) and other fine structure lines such as [SIV], [Nil], [NIII], 
and [OIII]. These line luminosities are diagnostic of key parameters such as the EUV luminosity 
and spectral shape, the X-ray luminosity and spectral shape, and the wind mass loss rate and shock 
speed. Our main results follows: 

1. The luminosity of fine structure lines (e.g., [Nell] and [Aril]) from the dominant ionization 
state of a species roughly scale with Lx and Lew- At very high Lx or Lew the lines saturate 
because the electron density in the emitting region exceeds the critical density of the line. [Aril] 
7.0 //m, which has not yet been observed, is predicted to be as strong as [Nell] 12.8 //m in the 
EUV layer. If the X-ray layer dominates and the X-ray spectrum is such that much of the X-ray 
luminosity is in the 1-3 keV band, the [Aril] line is predicted to be about 2.5 times weaker than 
the [Nell] line. Therefore, the observed [Nell] / [Aril] flux ratio may help determine the origin of 
these lines. Observations of [Nell] set upper limits for the EUV luminosity of the central star, 
^EUV ^ 10^^ EUV photons s~^ for most sources. 

2. Most of the fine structure emission in the EUV layer arises from 5-10 (M^/Mq) AU. Most 
of the fine structure emission from the X-ray layer is distributed more broadly in r from < 1 — 10 
AU for a solar mass star. 

3. If Lx ^ LeuVj there is about 2 times as much [Nell] emission arising from the X-ray layer 
as from the EUV layer, assuming our standard "soft" EUV (30,000 K blackbody) spectrum which 
produces the most [Nell] luminosity. 

4. A power-law EUV spectrum, L^, oc results in a [Nelll] line luminosity that is greater 
than the [Nell] line luminosity from the EUV layer, contrary to observations. If the EUV layer is 
responsible for the [Nell] emission, the EUV spectrum must be softer than a ~ 30, 000 — 40, 000 K 
blackbody spectrum between 15 eV and 40 eV. The X-ray layer, which has much higher abundances 
of atomic hydrogen, naturally gives [Nelll] line luminosities that are less than 0.1 of the [Nell] 
luminosities because of rapid charge exchange reactions of Ne"'"^ with H. 

5. Internal shocks in protostellar winds may be a viable explanation of the observed [Nell] in 
a number of sources, especially those with high M^, or its surrogate Mace- Confirmation of this 
origin requires high spatial ( ^ 1") and spectral ( ^ 10 km s~^) observations. The [Nell] from these 
regions, if they are nearby, may be extended (>1") and should produce broader (~ 100 km s~^ 
FWHM) profiles than the [Nell] from the EUV or X-ray layer, especially in face-on disks. 
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6. [01] 6300 A is weak [L^qi] ^ 10~^) from the EUV layer, the transition layer between the 
EUV layer and the X-ray layer, the X-ray layer if the spectrum is dominated by 1-2 keV photons, 
and likely also from the shear layer where the protostellar wind impacts the disk surface. A soft 
X-ray spectrum {L^ cx for 0.1 keV < hu < 2 keV) with considerable luminosity in 0.1 - 0.3 
keV photons produces a hotter and more ionized X-ray layer, and substantially more [01] 6300 
A luminosity because of the extreme temperature sensitivity of this line. Lx as high as 10~^ Lq 
with this spectrum results in L^qi] 10~^ Lq. The observed values of the low velocity component 
of [01] range from 10~^ to 10"'^ L0, with typical values ~ 10~^ L0. Therefore, soft X-rays are a 
plausible origin for the low velocity [01] component in many sources. 

7. We compared our models with a compilation of 54 sources of [Nell] emission from young 
low mass protostellar sources and with correlations of ^[Are//] with Lx and Mace- We note in point 
5 that internal shocks in winds may be a viable explanation for especially the sources with observed 
outflows or jets. There are also sources with low Mace where our "harder" X-ray spectrum, with 
most luminosity emerging at 1-2 keV, can explain the observed [Nell] emission. In some cases, the 
lines are resolved to be relatively narrow (10-60 km s^^), further indicating an X-ray layer origin 
and not a shock origin. However, there exist sources where neither wind shocks nor 1-2 keV X-rays 
carry sufficient energy to power the observed [Nell] line. These sources are likely candidates for 
[Nell] originating from the EUV layer or from an excess of soft (~ 0.1 — 0.3 keV) X-rays. If the 
spectrum in the EUV-soft X-ray wavelength region is a power law Li, oc i^^^, as Ercolano et al 
(2009) suggest, then the soft X-ray layer will dominate the production of [Nell], although the EUV 
layer may produce more [Nelll] than the X-ray layer. Whichever layer dominates, the [Nell] and 
[Nelll] luminosities directly provide a measure of the heretofore unobserved EUV or soft X-ray 
luminosities from the protostar or its immediate environs. 
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Table 1. IR Fine Structure Parameters for Species in 
Ionized Gas 



Transition 


C 
(er 
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gs) 


IT'ecr.s 

(cm~^) 
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Refs. 


[Aril] 7 fim 


1.8( 


-13) 


4.2(5) 


6.3( 


-6)^ 
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[Arlll] 9 


3.7( 


-14) 


1.2(6) 


6.3( 


-6)^ 


2 


[Nil] 122 


2.3( 


-14) 


1.6(3) 


9.1( 




3 


[NIII] 58 fim 


5.0( 


-13) 


1.2(3) 


9.1( 




4 


[Nell] 12.8 i^im 


2.2( 


-13) 


6.3(5) 


1.2( 


.4)b 


5 


[Nelll] 15.0 i^im 


3.7( 


-13) 


2.7(5) 


1.2( 


-4)b 


6 


[OIII] 52 


8.3( 


-13) 


4.6(3) 


3.2( 


-4)a 


3 


[SIII] 18 fim 


3.3( 


-13) 


1.5(4) 


7.6( 


-6)^ 


7 


[Sill] 33 fxm 


3.0( 


-13) 


4.1(3) 


7.6( 


-6)^ 


7 


[SIV] 10.4 fim 


9.5( 


-14) 


4.0(5) 


7.6( 


-6)^ 


8 



''^Abundances are from a. Savage & Sembach 1996, b. 
Grevesse &; Sauval 1998 and c. Asplund et al. 2005 
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Fig. 1. — The dependence of [Nell] 12.8 /xm and [Nelll] 15.5 luminosity with the EUV luminosity 
(top in erg s~^ and bottom in EUV photons s""*^) of the central star. The EUV spectrum is 
assumed to be a power law, Leuv{i^) oc u"^. This relatively hard EUV spectrum leads to high 
abundances of Ne^"*" in the EUV layer, and [Nelll] stronger than [Nell]. In Section 3 we explain the 
overall dependence of line luminosity proportional to EUV luminosity with saturation occurring 
at the higher luminosities as electron densities exceed the critical density of the [Nell] and [Nelll] 
transitions. 
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Fig. 2. — The dependence of [Nell] 12.8 fim and [Nelll] 15.5 luminosity with the EUV luminosity 
(top in erg and bottom in EUV photons s~^) of the central star. The EUV spectrum is assumed 
to be a blackbody with effective temperature T^jf = 30, 000 K. This relatively soft EUV spectrum 
leads to high abundances of Ne"*" in the EUV layer, and [Nell] significantly stronger than [Nelll]. 
In Section 3 we explain the overall dependence of line luminosity proportional to EUV luminosity. 
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Fig. 3. — The dependence of [Aril] 7 ^tm, [Arlll] 9 /.tm, [SIII] 19 fim, and [SIII] 33 fim line 
luminosities on the EUV luminosity (top in erg s~^ and bottom in EUV photons s~^) of the central 
star. The EUV spectrum is assumed to be a blackbody with effective temperature T^fj = 30, 000 K. 
Other lines from ionized species which require > 13.6 eV for their ionization are significantly weaker 
(see Table 1 and Eq. 8). We discuss analytic approximations for these predicted line luminosities 
in Section 3.1. 
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Fig. 4. — The EUV-produced line flux emergent from one side of the disk times 47rr^ is plotted 
against the radius of the disk. This luminosity is approximately the luminosity emerging from 
both sides of an annulus between 0.5r and 1.5r. The figure shows that most of the luminosity is 
generated at r ~ 10 AU. Shown are the results for a central star with Lew — Lx = 2 x 10'^'^ 
erg s~^. In photon units, ^euv = 10^^ s~^. The EUV spectrum is assumed to be a blackbody 
with effective temperature of 30,000 K (same case as Fig. 2). In addition, we have plotted (dotted 
line) 47rr^ times the emergent flux of [Nell] from the X-ray layer for our standard (" hard" ) X-ray 
spectrum. Substantial luminosity emerges from the region < 1 to 10 AU, and the overall [Nell] 
luminosity is ~ 2 times greater than the EUV layer in this case. 
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Fig. 5. — The gas temperature Tgas, the dust temperature T^ust, and the hydrogen nucleus density 
n are plotted versus the vertical distance z (top) from the midplane, or the hydrogen nucleus 
column N (bottom) from the surface (z = r). This vertical slice is for r = 10 AU. The central star 
X-ray luminosity and EUV luminosity and spectrum are the same as in Figure 4. Note that Tgas 
tracks Tdust to z ~ 2 AU, or ~ 10^^ cm~^. Higher in the disk, the gas is hotter than the dust. 
The EUV and X-ray layers are marked. Note that the ionization front is at zjp ~ 7.5 AU. Dust 
dominates the heating of the gas near the midplane. 
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Log N (cm-2) 

Fig. 6. — The fraction of neon which is in the singly ionized state, f{Ne'^), and the electron 
abundance relative to hydrogen nuclei, Xg, are plotted versus the hydrogen nucleus column N from 
the disk surface at r = 10 AU. The EUV layer only extends to about N ~ 10^^ cm~^ (see Figure 5) 
so that we highlight here the X-ray layer. The X-ray spectrum is our harder spectrum which peaks 
at ~ 2 keV. X-rays maintain the high f{Ne~^) throughout the region plotted. X-rays maintain a 
relatively high electron abundance to ~ 10^*^ cm~^. At higher columns, FUV photoionization of 
carbon as well as X-rays maintain Xg. 
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1028 1029 1030 1031 1032 

Lx (erg/s) 

Fig. 7. — The dependence of [Nell] 12.8 /im and [Nelll] 15.5 fim luminosity with the X-ray lu- 
minosity. The solid labelled lines are [Nell] and [Nelll] for our "harder" X-ray spectrum, where 
Ly (X V for 0.1 keV< hv < 2 keV. The dashed line is the [Nell] luminosity for our softer X-ray 
spectrum source, where Li, oc i^^^ for 0.1 keV < hu < 2 keV. Note the nearly linear dependence 
L[^e//] with Lx- Comparison with Fig. 2 shows that if Lx ~ Lew-, ^-nd assuming a soft EUV 
spectrum that produces the maximum amount of [Nell], then the [Nell] line luminosity is still 2 
times stronger from the X-ray layer as from the EUV layer. Also plotted are a recent compilation 
of [Nell] and Lx data (Giidel et al 2009). The region shaded with vertical dotted lines are sources 
with known outflows or jets. The region with horizontal solid lines are sources with undetected 
outflows. It appears that there are a substantial number of sources, especially the "outflow /jet" 
sources that are more luminous in [Nell] than the X-ray layer (or the EUV layer) could provide; 
internal shocks in the winds or jets are a possible explanation for these sources. 
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Fig. 8. — The dependence of [Nell] 12.8 /xm luminosity with the mass accretion rate onto the 
central star. The shaded regions are the same notation as in Figure 7; data from Giidel et al 
(2009). Section 3.4 in the text and Eq. (33) predict the [Nell] luminosity as a function of the wind 
or jet mass loss rate and the fraction fsh of the wind or jet that shocks at speeds greater than 
about 100 km s~^. The solid line in the figure assumes = O.lMacc and fsh = 1- The upper 
dashed line assumes = Mace and fsh = 1- The lower dashed line assumes that the product 
fsh^w is 10 times less than assumed in the solid line case. Shocks appear viable explanations for 
the origin of many of the [Nell] sources, especially those with observed outflows and jets (see also, 
Figure 7). 
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Fig. 9. — The predicted [01] 6300 A luminosity from the EUV layer is plotted for both a blackbody 
EUV spectrum {T(,ff = 30,000 K) or for a power law spectrum as a function of the EUV photon 
luminosity ^euv- The harder spectrum produces more [01] luminosity because more atomic O 
survives in the mostly ionized EUV layer (see text). Observed [01] luminosities are typically much 
higher than the ^ 10~^ Lq predicted from the EUV layer (see text). 
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Fig. 10. — The predicted [01] 6300 A luminosity from the X-ray layer is plotted versus the X-ray 
luminosity of the central star. The solid line is for our harder X-ray spectrum whereas the dashed 
line is for our softer X-ray spectrum (see text or caption to Fig. 7). Observed [01] luminosities 
are typically much higher than the < 10~^ Lq predicted from the X-ray layer produced by the 
harder spectrum. However, the softer X-ray spectrum produces [01] luminosities much more in 
accord with observations, because the X-ray layer is warmer and the line is extremely temperature 
sensitive (see text). 



